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13.  AOSTKACT  (MaoniMniOOwoffiK)  I 

Arclength  continuation  methods  were  incorporated  into  a  code  for  predicting  the  structure  of  sooting,  opposed-jet  flames. 
The  code  includes  complex  chemistry,  detailed  particle  dynamics,  particle  chemistry  and  radiation.  The  code  was  used  to 
predict  soot  production  over  a  wide  variation  in  strain  rates  for  both  ethylene/air  and  methane/air  diffusion  flames.  Predicted 
values  (both  peak  and  spatial  distributions)  agree  well  with  experimental  measurements  in  ethylene  flames.  Particle  size 
distributions  are  also  predicted  using  the  aerosol  equations  from  MAEROS,  but  no  data  is  available  for  comparison.  Also,  the 
soot  dynamical  equations  were  imbedded  into  a  separate  code  to  describe  soot  production  in  a  coflow,  laminar,  diffusion 
flame  which  includes  treatment  of  detailed,  gas  phase  chemistry.  Predictions  were  compared  to  measurements  made  in  a 
methane,  coflow  flame.  Reasonable  agreement  between  the  predictions  and  measurements  was  obtained,  although  a  factor  of 
three  underprediction  of  the  soot  volume  fractions  is  likely  due  to  uncertainties  in  inlet  conditions  and  an  inability  to  match 
closely  bulk  flame  parameters  such  as  temperature.  Predicted  peak  soot  production  occurred  around  1720K  and  particle! 
oxidation  was  dominated  by  superequilibrium  concentrations  of  hydroxyl  radicals.  Several  PAH-forming  sequences  were 
examined  and  compared  to  the  traditional  acetylene-addition  sequence.  A  sequence  involving  benzyl-propargyl  combination 
was  found  to  compete  with  the  traditional  mechanism  and  it  should  be  included  in  future  analyses.  The  algorithms  for 
treating  sectional  soot  dynamics  and  growth/oxidation  rates  were  modified  to  include  effects  at  high  pressure.  Continuum 
effects  and  limitations  to  gaseous  diffusion  were  included  in  the  opposed  jet  code.  Predicted  variations  in  soot  production  due 
to  pressure  changes  from  4  to  10  atmospheres  were  made  for  an  ethylene-air,  opposed  jet  flame. 
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I.  Objectives 

The  overall  objectives  of  this  work  are  to  (1)  refine  and  update  a  pyrolytic  model  which 
describes  PAH  formation;  (2)  modify  and  validate  an  opposed  jet  diffusion  flame  code 
which  includes  soot  formation  and  radiation  and  then  extend  model  predictions  over  a 
range  of  flame  parameters,  and  (3)  incorporate  the  soot  formation  and  radiation  models 
into  a  code  for  a  coflow  laminar,  diffusion  flame.  These  efforts,  particularly  the  latter 
two,  are  being  completed  with  assistance  firom  Prof.  M.  Smooke  of  Yale  University. 
Experiments  and  modeling  on  alkylated  aromatics  will  help  to  minimize  discrepancies 
between  known  PAH  formation  rates  and  those  required  to  explain  particle  inception 
rates  in  flames.  Validation  (and  modification  as  necessary)  of  the  existing  opposed  jet 
flame  code  wifii  data  from  ethylene  flames  will  increase  its  reliability  in  predicting 
sooting  conditions  with  other  hydrocarbon  fuels  and  inlet  conditions.  In  addition,  this 
effort  should  demonstrate  the  code’s  utility  prior  to  installing  soot  formation  and 
radiation  routines  into  codes  describing  other  flame  geometries  for  which  extensive  data 
bases  exist. 

II.  Status  of  Effort 

The  objectives  outlined  in  the  above  paragraph  have  been  met  during  the  course  of  this 
three  year  program.  Much  of  the  work  is  now  documented  in  publications  or  recently 
submitted  manuscripts;  including  (1)  Predictions  of  Soot  Dynamics  in  Opposed  Jet 
Diffusion  Flames  [1],  (2)  AppUcation  of  Continuation  Methods  to  Soot  Formation  in 
Counterflow  Diffusion  Flames  [2],  (see  Appendix  A),  (3)  Computational  and 
Experimental  Study  of  Soot  Formation  in  a  Coflow,  Laminar  Diffusion  Flame  [3],  (see 
Appendix  B)  and  by  several  presentations  at  technical  meetings,  (for  the  most  recent 
presentations  see  the  section  on  meeting  interactions). 

Most  of  the  research  described  in  [1]  was  a  result  of  work  completed  in  an  earlier  contract 
and  hence  will  not  be  treated  here.  The  results  from  [2],  [3]  and  a  recent  abstract  [4]  are 
the  main  accomplishments  of  this  contract.  A  brief  summary  of  conclusions  from  these 
studies  are  provided  in  Sections  III.  A.  and  B.,  while  tihe  details  are  presented  in  die 
Appendixes  of  this  report.  In  addition,  description  of  some  very  recent  work  on 
PAH/inception  modeling  and  the  development  of  algorithms  for  the  treatment  of  sectional 
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soot  dynamics  at  high  pressure,  not  yet  documented  elsewhere,  is  provided  in  Sections 
III.  C.  and  D. 

III.  Accomplishments/New  Findings 

III.  A.  Continuation  Model  for  Sooting,  Opposed  Jet  Flame 

Detailed  flame  models  that  include  complex  chemistry  coupled  with  soot  formation, 
growth,  oxidation,  and  particle  dynamics  equations  with  radiation  effects  have  been 
developed  for  both  laminar  premixed,  burner-stabilized  flames  and  opposed  jet  flames. 
These  models  have  been  based  upon  AFOSR  sponsored  work  described  previously  [1,5]. 
The  opposed  jet  flame  model  has  been  extended  by  incorporating  an  arclength 
continuation  method  into  the  code  such  that  flame  solutions  can  be  readily  obtained  over 
a  wide  range  of  strain  rates.  This  recently  documented  procedure  [2]  dramatically 
decreases  the  effort  required  to  obtain  solutions  over  a  range  of  conditions.  Solutions 
were  obtained  from  low  strain  to  extinction  for  both  ethylene  and  methane  flames.  Peak 
soot  volume  fractions  were  found  to  decrease  gradually  (logarithmic  scale)  from  low  to 
high  strain  by  about  two  orders  of  magnitude  for  the  ethylene  flames.  The  decay  was 
principally  attributed  to  the  decrease  in  residence  time  in  the  high  strain  flames,  although 
decreasing  benzene  concentrations  and  flame  temperatures  contribute  to  the  decrease.  The 
concentration  of  acetylene,  the  principal  species  contributing  to  surface  growth,  did  not 
vary  by  more  than  a  factor  of  two.  Predicted  soot  volume  fractions  agreed  well  with 
experimental  measurements  of  peak  values  and  the  spatial  dependence  of  soot  volume 
fractions  for  the  ethylene  flames.  There  is  yet  no  data  available  to  compare  the  predictions 
of  particulate  size  and  how  the  distribution  varies  with  location  in  the  opposed  jet  flame. 

III.  B.  Soot  Formation  in  a  Coflow,  Laminar,  Diffusion  Flame 

In  a  second  study  [3],  the  soot  dynamical  equations  have  been  coupled  with  the  equations 
describing  an  axisymmetric,  coflow,  laminar  diffusion  flame.  Predictions  have  been  made 
for  a  methane  jet  flame  and  have  been  compared  to  recent  experimental  measurements. 
Effects  of  thermal  radiation  and  particle  scrubbing  of  the  gas  phase  species  (for  inception, 
growth,  and  oxidation)  are  included.  General  features  of  the  flame  are  reasonably  well 
described,  although  peak  soot  volume  fractions  are  underpredicted  by  about  a  factor  of 
three.  We  believe  ^s  underprediction  is  principally  due  to  the  inability  to  predict 
accurately  the  bulk  flame  properties,  such  as  low  predictions  of  centerline  temperatures 
and  a  slight  overprediction  of  the  flame  height.  Possible  sources  of  these  uncertainties 
include  radiative  heat  transfer  to  the  uiueacted  fuel,  burner  lip  effects,  and  perhaps 
insufficient,  low  temperature,  fuel  rich  chemical  kinetics  for  the  methane  fuel.  (As 
described  in  Section  III.C.,  there  may  also  be  uncertainties  in  the  calculated  inception 
rates.)  Predicted  distributions  of  particle  sizes  indicate  the  existence  of  relatively  small 
particles  along  the  centerline  of  the  burner  and  small  concentrations  of  large  particles 
which  escape  the  flame  front,  unoxidized.  Oxidation  of  particulates  is  dominated  by 
hydroxyl  radicals  which  attain  concentrations  about  ten  times  above  equilibrium  near  the 
flame  front.  Gas  cooling  effects,  due  to  radiative  losses,  are  shown  to  have  a  very 
significant  effect  (factor  of  three  decrease)  on  predicted  soot  volume  fractions,  for  this 
lightly  sooting,  methane  flame. 
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The  solution  file  for  this  coflow,  methane  flame  has  been  recently  interrogated  [4]  (see 
Appendix  C)  to  extract  information  about  the  soot  formation  process  in  the  flame.  Peak 
surface  growth  rates  were  found  to  be  on  the  rich  side  of  the  flame  front  at  temperatures 
of  1720-1760K,  which  is  in  the  vicinity  of  the  peak  of  the  traditional  bell-shaped  curve 
for  soot  formation.  Local  equivalence  ratios  under  these  conditions  are  about  1.8,  tending 
towards  less  rich  conditions  higher  in  the  flame. 


III.  C.  PAH  Formation  and  Modeling 

Previously,  we  have  argued  that  existing  detailed  inception  models  imbedded  in  soot 
formation  codes  generally  underpredict  the  formation  rates  of  PAH  molecules  in  flame. 
Hence  it  is  likely  that  predicted  inception  rates,  based  upon  traditional  routes  for  PAH 
production,  are  low.  While  there  are  some  arguments  to  the  contrary  [6],  we  have  yet  to 
adopt  the  higher  reaction  rates  assumed  in  the  alternate  analysis.  Instead,  we  believe  that 
there  are  mechanisms  contributing  to  PAH  production  in  addition  to  the  alternate 
acetylene  addition  steps  (and  occasional  benzene/phenyl  addition)  and  cyclization 
processes.  Based  on  the  recent  work  on  the  toluene  pyrolysis  mechanism,  with  and 
without  additives  [7,8],  we  have  just  considered  relative  inception  rates  based  on  the 
reaction  set  shown  in  Table  I. 


Table  I. 


QHs 

+ 

QHe 

+ 

C^Ks 

QHs 

+ 

CA 

H- 

C,H, 

+ 

C14H14 

+ 

Alternate  Mechanism  for  PAH  Production  via  Toluene/Benzyi 


CH3 

<=> 

QHs 

CH3 

C3H8 

+ 

H 

H 

CA 

+ 

H, 

CH3 

CA 

+ 

CH4 

C,H« 

CA 

+ 

H 

C3H3 

c,A 

+ 

2H 

QH, 

<=> 

Ci4H,4 

H 

=> 

C14H10 

+ 

2H2 

-f 

H 

log 

log 

(Afor) 

F 

•“^forward 

(ArJ 

F 

■^reverse 

12.92 

-800 

16.0 

97000 

12.3 

17000 

13.08 

5100 

14.1 

8400 

13.5 

23700 

12.2 

11100 

13.08 

27000 

15.45 

88900 

14.24 

-200 

11.7 

5000 

12.7 

-500 

14.93 

60000 

13.7 

5000 

*  A  -  cc,  moles,  sec  units;  E  -  calories/mole 


Forward  and  reverse  rate  constants  in  this  table  are  taken  from  the  previous  analysis  on 
toluene  pyrolysis  [7,9]  but  with  an  updated  rate  constant  [8]  for  the  propargyl-benzyl 
recombination  reaction  that  forms  naphthalene  (CioH*). 


An  analysis  based  on  the  steady-state  assumptions  for  toluene,  bibenzyl  and  the  benzyl 
radical  was  performed  to  estimate  a  rate  of  production  of  naphthalene  from  this  reaction 
sequence.  The  rates  of  production  of  naphthalene  via  this  mechanism  were  compared  to 
the  rates  of  production  of  CjoH,  via  the  more  traditional  acetylene  addition  mechanism. 
Using  the  local  solutions  of  species  concentrations  and  temperature  from  the  coflow. 
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methane  flame,  the  rates  from  this  ‘new’  mechanism  were  found  to  be  at  least  comparable 
to  those  for  the  traditional  route  in  the  flame  front  and  generally  nearly  an  order  of 
magnitude  higher  near  the  centerline  of  the  flame.  These  results  qualitatively  support  the 
experimental  result  obtained  by  [10]  that  addition  of  small  concentrations  of  toluene  in 
the  fuel  increases  the  production  of  naphthalene  in  the  coflow,  methane  flames. 

In  addition,  we  have  examined  the  potential  for  recombination/rearrangement  of 
cyclopentadienyl  radicals  as  a  source  of  naphthalene  in  these  flames.  The  kinetic 
calculations  indicate  that  this  route  should  dominate  the  naphthalene  production; 
however,  we  do  not  yet  endorse  this  route  as  the  experimental  evidence  for  this  process  is 
less  clear  and  it  is  likely  that  we  have  significantly  overestimated  this  net  rate  due  to  the 
lack  of  including  scrubbing  effects  and  the  skeletal  mechanism  utilized  for  Cj-chemistry. 

Nevertheless,  these  results  strongly  support  the  necessity  of  including  alternate  PAH 
formation  routes  in  the  modeling  of  inception  processes.  Such  sequences  have  not  yet 
been  included  in  our  flame  modeling  but  are  anticipated  to  be  installed  shortly.  The  role 
of  these  alternate  mechanisms  must  be  considered  more  completely  in  subsequent 
modeling  of  soot  formation  in  flames. 

III.  D.  Treatment  of  Sectional  Soot  Dynamics  at  High  Pressure 
The  sectional  soot  growth  program  has  been  modified  to  be  valid  for  high-pressure 
calculations.  The  first  versions  had  been  written  for  the  low-pressure,  free-molecule 
regime  in  which  mean  free  path  is  l^ger  than  characteristic  particle  sizes  (pressures  up  to 
approximately  1-2  atmospheres).  Coalescence  rates  in  the  free-molecule  regime  are  equal 
to  sticking  probabilities  times  hard-sphere  collision  frequencies,  surface  growth  and 
oxidation  growth  coefficients  have  a  simple,  semi-analytic  form,  and  thermophoretic 
velocity  is  independent  of  particle  size  [1].  As  pressure  increases,  this  assumption  of 
mean  free  path  larger  than  particle  sizes  becomes  invalid,  first  for  the  largest  particles, 
and  the  so-called  continuum  regime  is  entered.  The  intermediate  situation  between  the 
strict  free-molecule  and  continuum  cases  is  called  the  transition  regime.  Corrections  to 
the  rates  of  coalescence,  surface  growth  and  oxidation,  and  thermophoresis  are  then 
required. 

The  Sectional  coalescence  coefficients  were  modified  for  continuum  effects  using 
algorithms  already  contained  in  the  MAEROS  program  [11]  from  which  the  soot  program 
originated.  The  contributions  of  Brownian  motion  to  the  coalescence  rates  in  the 
transition  and  continuum  regime  are  treated  in  the  manner  of  Fuchs  [12].  Gravitational 
and  turbulence  contributions  were  not  included.  As  continuum  effects  become  important, 
the  thermophoretic  velocities  of  the  largest  particles  start  to  decrease,  so  that  in  general  a 
variation  of  velocity  with  size  class  must  be  considered.  At  very  high  pressure, 
thermophoretic  velocity  again  becomes  independent  of  size  class.  The  transition  and 
continuum  corrections  were  again  calculated  using  algorithms  embedded  in  the 
MAEROS  program,  which  were  in  turn  based  on  the  treatment  of  Waldmann  and  Schmitt 
[13]. 
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Continuum  effects  corrections  for  surface  growth  and  oxidation  are  more  complicated. 
The  net  rate  of  growth  or  oxidation  at  a  particle  surface  has  contributions  both  from 
diffusion  of  gaseous  reactants  to  the  surface,  and  the  intrinsic  chemical  process  at  the 
surface.  Measured  rates  are  true  surface  reaction  rates  only  if  they  are  small  relative  to 
the  rate  for  diffusion  of  reactants  to  the  surface.  In  the  free-molecule  limit,  the  reactant 
flux  at  the  surface  is  given  by  gas  kinetic  impingement  rates,  and  reactant  concentrations 
at  the  surface  are  substantially  equal  to  the  asymptotic  values  far  from  the  surface.  As 
the  pressure  increases  to  the  continuum  limit,  diffusion  theory  is  used  to  calculate  the 
reactant  flux  [14].  The  transition  regime  diffusional  reactant  flux  is  calculated  using  an 
approximate  interpolation  formula  proposed  by  Fuchs  and  Sutugin  [14,15].  With  the 
diffusional  flux  at  the  particle  surface  established,  the  assumption  of  steady  state  at  the 
surface  leads  to  resistance  rule  expressions  for  the  net  rate  of  carbon  atom  addition  due  to 
growth,  and  carbon  atom  removal  due  to  oxidation  by  oxygen  and  OH.  This  expression 
has  the  proper  limiting  forms.  In  the  limit  that  the  diffusional  flux  is  much  larger  than  the 
intrinsic  surface  processes,  the  total  per  particle  growth/oxidation  rates  will  always  be 
proportional  to  particle  surface  area,  or  the  mass  of  the  particle  raised  to  the  two-thirds 
power.  When  the  intrinsic  rates  dominate  in  the  free-molecule  regime,  the  proportionality 
to  surface  area  will  again  apply.  When  the  intrinsic  rates  dominate  in  the  continuum 
limit,  the  total  per  particle  rate  will  be  proportional  to  particle  diameter,  or  to  the  square 
root  of  particle  surface  area.  In  the  high-pressure  soot  program,  the  generalized  rate 
expression  is  employed  in  the  numerical  integration  over  the  assumed,  intra-sectional 
mass  distribution  function  which  is  used  to  calculate  the  sectional  coefficients  [1,11]. 
This  is  done  individually  for  the  acetylene  addition  growth  process,  oxidation  by  oxygen, 
and  oxidation  by  OH.  One  assumption  of  the  present  program  is  that  for  simplicity  the 
algorithm  has  been  set  up  assuming  intrinsic  rates  first  order  in  the  applicable  reactant 
(acetylene,  oxygen,  OH).  Particle  diffusional  transport  makes  such  a  small  contribution 
in  most  simulations  that  for  simplicity  the  particle  diffusional  coefficients  have  been  set 
equal  to  zero  in  the  high  pressure  program.  The  program  calculates  all  coefficients  at  a 
given  pressure  for  a  range  of  temperafiires,  and  then  uses  interpolation  to  extract 
coefficients  at  intermediate  temperatures. 

Figure  1  shows  an  example  calculation  for  a  plug  flow,  ethylene-air  opposed  jet  at  a  str^ 
rate  of  80  sec'*.  The  sensitivity  to  pressure  is  highest  at  the  lowest  pressures,  becoming 
less  sensitive  at  higher  pressures  as  the  importance  of  particle  scrubbing  of  gaseous 
growth  species  becomes  important.  The  ability  to  account  for  continuum  effects  on  the 
particle  growth  and  transport  coefficients  may  now  be  on  a  sounder  basis  in  this  code 
than  is  the  ability  to  model  pressure  effects  on  the  gas  phase  chemistry. 
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Figure  1.  Predicted  Sensitivity  of  Soot  Volume  Fraction 
to  Pressure  in  Ethylene-Air  Opposed  Jet  Flame 
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Abstract 


An  arclength  continuation  method  has  been  incorporated  into  a  coupled  gas/particle  op¬ 
posed  jet  soot  prediction  code  and  used  to  examine  the  soot  production  in  methane-air 
and  ethylene-air  diffusion  flames.  The  model  includes  a  detailed  transport,  flmte  rate 
chemistry  gas  phase  system  coupled  with  aerosol  dynamical  equations  in  the  sectional 
representation.  The  problem  formulation  includes  detailed  submodels  for  the  inception, 
growth,  oxidation,  and  dynamics  of  soot  particulates.  The  continuation  feature  produces 
a  dramatic  reduction  in  the  effort  required  to  carry  out  extensive  studies  of  the  sensitivity 
of  soot  growth  to  flame  parameters  sudi  as  the  strain  rate.  Computations  for  both  atmo¬ 
spheric  pressure  methane-air  and  ethylene-air  counterflow  diffusion  flames  are  presented. 
Effects  of  radiation  and  variation  in  important  soot  growth  species  are  examined.  Pull 
coupling  of  radiation  and  scrubbing  within  the  flame  and  soot  equations  is  demonstrated 
to  be  dramatically  important,  especially  at  low  strain  rates.  Residence  time  at  elevated 
temperature  was  found  to  be  a  dominant  parameter  controlling  soot  production. 
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APPLICATION  OF  CONTINUATION  METHODS 
TO  SOOT  FORMATION  IN  COUNTERFLOW  DIFFUSION  FLAMES 

Introduction 

Numerical  models  for  analyzing  the  structure  of  opposed  jet  diffusion  flames  (See  Fig¬ 
ure  1)  with  detailed  chemistry  have  been  refined  over  the  past  ten  years  [1-5].  Current 
modeling  capabilities  for  opposed-jet,  non-sooting  flames  are  now  excellent,  with  most  un¬ 
certainties  centered  on  inlet  conditions  and  on  the  chemical  kinetic  sets  used  to  model  the 
data.  Comparisons  between  experimental  and  computational  results  have  been  made  with 
the  goal  of  helping  to  interpret  flame  structure,  extinction  limits  and  soot  formation  phe¬ 
nomena.  In  particular,  there  have  been  studies  on  soot  inception  limits  in  air  [6-7],  effects 
of  oxygen  in  the  fuel  [8],  effects  due  to  strain  rate  variations  [9],  and  on  the  relative  impor¬ 
tance  of  temperature  and  fuel  dilution  on  soot  formation  processes  [10].  Typically,  only 
skeletal  models  integrating  the  soot  formation  processes  with  opposed-jet  flame  computa¬ 
tions  have  been  available  to  assist  in  the  interpretation  of  the  experimental  data.  Recently, 
a  comprehensive  particle  kinetics  model  was  constructed  [11]  which  provided  qualitative 
predictions  comparable  to  existing  data  for  soot  formation  in  an  opposed-jet,  methane-air 
diffusion  flame  [12].  Particle  sectional  equations  with  submodels  for  particle  inception, 
growth,  coagulation,  and  oxidation  [13]  were  coupled  to  the  gas  phase  conservation  equa¬ 
tions  for  an  opposed-jet  diffusion  flame.  Complex  chemistry  and  detailed  transport  along 
with  particulate  and  gas  radiation  were  included.  The  model  has  the  capability  to  provide 
detailed  spatially-dependent  information  on  volume  fractions,  particle  number  density,  size 
and  size  distribution.  Codes  of  this  type  can  be  used  to  validate,  over  required  ranges  of 
temperatures,  pressures  and  fuel  types,  various  proposed  (semi-empirical)  mechanisms  for 
modeling  soot  formation  and  to  help  determine  simplifying  assumptions  which  can  be 
made  in  flmd-dynamically  more  complex  models.  The  model  developed  in  [11]  has  been 
exercised  recently  in  an  attempt  to  compare  several  soot  formation  mechanisms  proposed 
and  to  re-evaluate  some  of  the  temperature  and  dilution  effects  described  in  the  literature 
[14]. 

Solutions  of  strained,  nonpremixed  laminar  flames  have  been  used  in  the  modeling  of 
turbulent  reacting  flows  in  the  laminar  flamelet  regime  [15].  Such  apphcations  require 
solutions  over  a  large  range  of  strain  rates,  which  generate  the  functional  relationships 
between  soot  volume  fraction  and  mixture  fraction/ scalar  dissipation.  Mohual  generation 
of  these  solutions  is  highly  labor  intensive.  To  resolve  this  problem,  the  model  derived  in 
[11]  has  been  modified  to  include  continuation.  By  incorporating  arclength  continuation 
concepts  [16,17],  one  can  now  perform  accurate,  highly  resolved,  detailed  parameter  studies 
in  which  the  strain  rate,  pressure,  inlet  temperatures,  dilution  factor  and  (for  premixed 
cases)  the  equivalence  ratio  are  varied.  The  model  and  solution  algorithm  are  described 
in  this  manuscript.  Predictions  for  soot  volume  fractions  as  a  function  of  strain  rate  have 


2 


been  made  and  reported  herein  for  both  ethylene-air  and  methane-arr  diffusion  flames  at 
atmospheric  pressure  with  the  reactants  at  room  temperature.  Effects  of  radiation  are 
examined  and  soot  particle  dynamics  are  analyzed. 

Soot  Modeling 

The  soot  kinetics  continuation  analysis  is  based  on  a  single  solution  opposed  jet  particle 
growth  program  that  employs  the  sectional  representation  for  spheroid  growth  [11].  The 
picture  of  soot  kinetics  is  that  of  coalescing,  solid  carbon  spheroids  undergoing  surface 
growth  in  the  free  molecule  Umit.  The  particle  mass  range  of  interest  is  divided  into 
sections  [18],  and  a  dynamical  equation  written  for  each  section  including  coalescence, 
surface  growth,  and  oxidation.  For  the  first  bin,  an  inception  source  term  is  included.  The 
incorporation  of  these  dynamical  equations  into  a  transport/conservation  equation  for 
each  section  includes  both  thermophoresis  and  an  effective  bin  diffusion  rate.  In  the  gas 
species  conservation  equations,  provision  is  made  for  complex  chemistry;  additional  source 
terms  are  included  to  account  for  scrubbing  or  generation  of  gaseous  species  arising  from 
the  particle  growth  and  oxidation  processes.  The  gas  and  soot  are  additionally  coupled 
through  non-adiabatic  radiative  loss  from  both  the  gas  and  the  soot  in  the  optically  thin 
approximation.  For  one  atmosphere  calculations,  the  sectional  coefficients  for  transport, 
coalescence,  and  surface  processes  are  calculated  in  the  free  molecule  limit  [18],  with  surface 
growth  and  oxidation  proportional  to  particle  surface  area.  In  the  free  molecule  hmit,  the 
particle  mass  and  medium  dependencies  of  the  sectional  coefficients  are  factored  in  such  a 
way  that  they  only  have  to  be  evaluated  once.  Oxidation  of  soot  by  O2  and  OH  is  treated 
as  described  in  [11].  The  inception  model  employed  here  is  based  on  an  estimate  of  the 
formation  rate  of  two-and  three-ringed  aromatic  species,  and  is  a  function  of  local  acetylene, 
benzene,  phenyl,  and  molecular  hydrogen  concentrations.  Assuming  steady-state  values  of 
intermediates  and  that  »  C2H2,  the  rate  of  production  of  the  polyaromatic  species 
can  be  estimated  to  be  [11] 


4C10H7]  ^  ^Qii.88g(-4378/T)I^^jC'6if5]  cc/mole/sec,  (1) 

dt  [Hi] 

and 

d[CuHio]  ^  ^Qn.50^(-e390/T)\22MA[C^He][CeH^]  cc/mole/sec,  (2) 

dt  [Hi] 

where  the  gas  phase  concentrations  and  temperatures  are  evaluated  at  local  conditions. 
With  the  further  assumptions  that  inception  is  limited  by  the  formation  of  polyaromatics 
and  oxidation/decomposition  of  such  species  can  be  neglected,  the  inception  rate.  Si,  in 
grams/cc/sec  was  initially  assigned  to 


3 


(3) 


Si  =  127  X 


d[CroHr] 

dt 


+  178X 


d[C^M 

dt 


where  the  constants  (molecular  weights)  are  provided  to  convert  from  molar  to  mass  units. 
This  inception  rate  is  included  as  a  source  term  in  the  dynamical  equations  for  the  first 
size  class,  whose  lower  mass  boundary  is  set  at  127  amu.  The  upper  mass  boundary  of 
the  highest  class  is  set  sufficiently  large  to  ensure  that  the  solution  is  independent  of  this 
maximum  massj  for  these  calculations  the  maximum  value  is  approximately  10*^^  amu. 
There  is  no  significant  density  of  particles  in  the  highest  classes. 

The  surface  growth  rate  is  nominally  that  of  Harris  and  Weiner  [19],  who  report,  from 
premixed  flame  investigations,  a  rate  proportional  to  acetylene  concentration;  their  rate  is 
normalized  in  this  application  to  1650K,  with  a  31.8  kcal  activation  energy  as  suggested  by 
Hura  and  Glassman  [8].  However,  there  is  evidence  from  ethylene  diffusion  flames  that  the 
effective  surface  growth  rate  is  somewhat  larger  than  that  reported  by  Harris  and  Weiner, 
who  used  optical  techniques  and  the  assumption  of  isolated  soot  spheroids  (Axelbaum, 
Flower,  and  Law  [20];  Kennedy,  Kollmann,  and  Chen  [21];  Sunderland  and  Faeth  [22]). 
Puri,  Richardson,  Santoro,  and  Dobbins  [23]  have  suggested  that  the  assumption  of  spher¬ 
ical  particle  shapes  in  analysis  of  light  scattering  data  can  lead  to  error,  and  that  a  more 
appropriate  interpretation  of  the  data  in  terms  of  an  aggregate  particle  morphology  can 
be  made,  leading  to  different  total  soot  surface  areas  and  specific  growth  rates.  As  will 
be  seen,  it  has  been  found  on  an  empirical  basis  that  it  is  necessary  to  enhance  the  HW 
growth  rate  by  about  a  factor  of  two. 

The  main  soot  parameter  of  interest  in  this  paper  is  the  soot  volume  fraction,  for  which 
convergence  is  generally  obtained  with  a  smaller  number  of  sections  than  is  the  case  for 
particle  size  and  number  density  [11].  The  munber  of  sections  required  for  convergence 
must  be  examined  in  each  case,  although  evidence  points  to  the  relative  magnitudes  of 
surface  growth  and  inception  as  the  important  parameter.  When  the  contribution  of 
inception  to  total  volume  fraction  is  significant,  although  not  necessarily  dominating,  the 
number  of  sections  required  can  be  very  small  [11],  with  larger  numbers  necessary  when 
the  contribution  of  surface  growth  to  volume  fraction  becomes  very  much  larger  than  that 
from  inception.  A  total  of  twenty  soot  classes  were  used  in  the  calculations  reported  in  this 
paper.  This  level  insured  that  predicted  soot  volume  fractions  were  essentially  independent 
of  the  number  of  sections. 

The  chemical  mechanism  used  for  modeling  the  methane  flames  is  a  modified  version 
of  GRIMech  2.11  [24]  with  1)  the  pressure-dependent  reactions  evaluated  at  atmospheric 
pressure,  2)  the  NOx  chemistry  removed  and  3)  the  addition  of  the  formation  and  oxidation 
steps  of  benzene  and  related  species  [11].  It  contains  38  species  and  192  individual  reaction 
steps.  The  ethylene  species  mechanism  employed  has  45  species  and  233  reactions.  This 


4 


set,  derived  from  GRIMech  1.2,  has  been  developed  at  UTRC  based  upon  comparisons 
to  experimental  data  from  perfectly  stirred  and  flow  reactors  as  well  as  ignition  delay 

measurements. 

As  an  example  of  how  the  model  predictions  agree  with  experimental  data,  the  single 
solution  code  has  been  used  to  simulate  the  ethylene  opposed  jet  experiments  of  Wang  et 
aJ.  [25].  Although  the  experiment  is  nominally  described  by  plug  flow,  Wang,  et.  al.  report 
that  the  real  flow  is  a  mixture  of  plug  and  potential  solutions,  and  describe  a  process  by 
which  they  adjust  boundary  conditions  to  match  their  experimental  velocity  profiles  in  the 
flame  region  using  a  potential  flow  simulation.  The  simulations  described  here  employed 
a  different  procedure.  Plug  flow  is  assumed,  but  the  velocity  on  the  fuel  side  is  adjusted 
to  bring  the  theoretical  velocity  profile  into  agreement  with  experiments.  The  oxidizer 
velocity  is  assigned  the  value  determined  from  the  stated  oxidizer  flow  rate.  Figure  2 
shows  that  this  procedure  gives  excellent  agreement  with  data  for  the  highest  and  lowest 
strain  rates  reported  by  Wang,  et.  al..  With  such  fits,  it  can  be  assumed  that  the  particle 
residence  times,  important  for  soot  growth  simvJations,  are  accurate  representations  of  the 
experiments.  The  velocity  gradient  shown  is  the  difference  of  the  fuel  and  oxidizer  velocities 
divided  by  the  distance  between  the  fuel  and  oxidizer  jets  required  to  fit  the  velocity  data. 
Figure  3  shows  the  comparisons  of  predicted  and  measured  soot  volume  fractions  for  an 
assumed  factor  of  two  enhancement  of  the  HW  surface  growth  rate.  Deviations  from  the 
experimental  results  are  well  within  experimental  uncertainty. 

The  surface  growth/inception  models  used  here  to  illustrate  the  soot  continuation  anal¬ 
ysis  agree  well  with  available  ethylene  opposed  jet  data.  Nevertheless,  it  is  not  suggested 
that  it  represents  a  universal  solution  to  the  problem  of  soot  formation  because  compar¬ 
isons  with  soot  data  in  other  fuels  have  not  been  carried  out.  The  possibility  of  aging 
of  the  surface  growth  rate  has  not  been  accounted  for  here,  nor  are  the  results  with  fi¬ 
nite  coalescence  appropriate  for  aggregate  formation.  Fundamental  questions  relating  to 
inception,  temperature  dependence  of  surface  growth,  and  precise  dependence  of  surface 
growth  on  gas  phase  species  concentrations  remain.  The  need  to  consider  gaseous  PAH  as 
major  surface  growth  species  in  addition  to  acetylene  has  also  been  suggested  by  Benish, 
Lafleur,  Taghizadeh  and  Howard  [26].  To  some  extent,  such  processes  are  included  m  the 
formulation  of  our  model,  since  PAHs  are  implicitly  treated  as  soot  particles  and  hence 
are  included  in  the  coalescence  and  siufface  growth  equations.  Whether  soot  growth  rates 
derived  from  premixed  and  opposed  jet  studies  are  self-consistent  is  still  an  open  issue,  as 
well. 

Problem  Formulation 


Conservation  Equations 

To  model  the  combined  gas  phase/soot-laden  counterflow  diffusion  flame  (Figure  1), 
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we  utilize  the  elliptic  form  of  the  two-dimensional  conservation  equations  in  cylindrical 
coordinates  for  the  gas  phase  together  with  corresponding  sectional  species  equations  (as¬ 
suming  no  slip  between  the  gas  and  particle  phases).  The  complexity  of  the  coupled  sys¬ 
tem  can  be  reduced  [27]  by  seeking  a  similarity  solution  of  the  form  u  =  rU(z),  v  =  v(z), 
Vk  =  Vkiz),  k  =  1,2,...K,K  +  +  M,T  =  T(z),  where  r  and  z  denote  the 

independent  radial  and  axial  spatial  coordinates,  respectively;  T,  the  temperature;  K, 
the  total  number  of  gas  phase  species  and  Af,  the  number  of  soot  (size)  sections;  Yt, 
the  mass  fraction  of  the  k*'‘  gas  phase  species  (k  <  K)  or  the  (k  —  K)*’'  soot  size  class 
(k  >  K')  and  u  and  t;  the  radial  and  axial  components  of  the  velocity,  respectively.  If  we 
substitute  the  expressions  for  the  two  velocities  into  the  continuity  and  the  momentum 
equations,  we  find  that  the  reduced  pressure  gradient  in  the  radial  direction  is  constant, 
i.e.,  Ifr  ■  dPidr  =  J,  where  P  denotes  the  pressure.  If  we  use  this  result,  together  with  the 
similarity  assumptions,  we  find  that  the  combined  gas  and  particle  two-dimensional  set  of 
governing  equations  can  be  reduced  to  a  nonlinear  two-point  boundary  value  problem  in 
the  axial  direction.  The  governing  equations  for  mass,  momentum,  chemical  species,  soot 
size  classes  and  energy  can  be  written  in  the  form 

^+2pU  =  0,  (4) 

dz 


d  (  dV\  ,.dU  ,  „ 

-^(.pYM.)  -  +  (*|  +  WDIV.  =  0,  h  =  1,2,...,  K, 

J  jyr 

+  Vn  J)  -  v-^  +  Q*  =  0,  k  =  K+YK  +  2,...,K  +  M, 

d  {,dT\  ^.dT  ^  dT  ^  ,  dT 


(5) 

(6) 
(7) 


K+M 

E 


(wl  +  wl)Wkhk 


+  ^  =  o, 

dz 


(8) 


and 


PW 


(9) 


where  V  =  pv,  and,  in  addition  to  the  variables  already  defined,  p  denotes  the  mass  density; 
Wjt,  the  molecular  weight  of  the  k*^  species;  W,  the  mean  molecular  weight  of  the  mixture; 
i?g,  the  universal  gas  constant;  A,  the  thermal  conductivity  of  the  mixture;  Cp,  the  constant 
pressure  heat  capacity  of  the  mixture;  Cpk,  the  constant  pressure  heat  capacity  of  the  k 
species;  twf,  the  molar  rate  of  production  of  the  fc*'"  species  per  unit  volume  due  to  gas 
phase  chemistry,  wl,  the  molar  rate  of  production  of  the  species  per  unit  volume  due 
to  scrubbing/replenishment  by  the  soot  growth/oxidation  processes;  Qk  the  rate  of  change 
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of  section  k  due  to  aerosol  processes;  ht,  the  specific  enthalpy  of  the  A:*'*  species;  /x  the 
viscosity  of  the  mixture;  Vfc.,  the  diffusion  velocity  of  the  A:*'-  species  in  the  direction;  ur 
and  Vd^  the  thermal  and  mass  diffusion  velocities  of  the  A:*'‘  soot  size  class  and  dq^J dz  the 
divergence  of  the  net  radiative  flux  for  gas  bands  and  soot  in  the  optically  thin  limit.  To 
complete  the  specification  of  the  problem,  boundary  conditions  must  be  imposed  on  both 
the  fuel  and  oxidizer  jets.  (For  more  details  the  reader  is  referred  to  [11]). 


It  is  worthwhile  to  note  that  in  the  above  formulation  the  strain  rate  does  not  appear 
explicitly  in  the  governing  equations.  The  magnitude  of  the  reduced  pressure  gradient,  J, 
is  related  to  the  “stretch”  in  the  flame  due  to  the  imposed  flow.  As  the  magnitude  of  J 
increases,  so  does  the  strain  rate.  The  strain  rate  can  be  computed  from  the  axial  velocity 
gradient.  However,  since  this  quantity  varies  throughout  the  flame,  one  must  decide  on 
a  specific  location  (for  example,  the  maximum  value  of  the  axial  velocity  gradient  on  the 
oxidizer  or  the  fuel  side  have  been  used)  for  its  evaluation.  Alternatively,  one  could  employ 
the  asymptotic  formula  derived  in  [28].  In  this  paper  we  utilize  an  “average”  strain  rate 


where  Vj  and  are  the  fuel  and  oxidizer  velocities,  respectively,  and  L  is  the  length  of  the 
computational  domain.  We  point  out  that  for  ethylene-air  diffusion  flames,  the  expression 
in  (10)  is  essentially  equal  to  one  half  the  value  defined  in  [28].  A  major  advantage  in  using 
(10)  is  that  it  can  easily  be  obtained  from  experimental  flow  rates  and  does  not  require 


velodty  measurements  inside  the  domain. 


Method  of  Solution 

Our  goal  is  to  be  able  to  study  the  variation  of  soot  volume  fraction  as  a  function  of  the 
strain  rate.  As  the  strain  rate  varies  from  low  to  high  values,  the  flame  nears  extinction. 
During  this  process  the  flame  can  adjust  its  position  in  the  computational  domain.  As 
the  extinction  point  is  reached,  the  Jacobian  matrix  becomes  singular.  This  precludes 
the  use  of  conventional  solution  algorithms.  As  a  result,  to  be  able  to  move  smoothly 
from  low  through  extinction  values  of  the  strain  rate,  we  employ  an  arclength  continuation 
solution  algorithm  with  Newton-Uke  iterations  and  global  adaptive  gridding  [17].  The  set 
of  equations  (4-9)  reduces  to  a  system  of  the  form 

.F(Af,J)  =  0,  ,  (11) 

where  X  =  Yk+m,  T)  is  the  solution  vector.  The  solutions  (X ,  J)  in  (11)  form 

a  one-dimensional  manifold  which,  as  a  result  of  the  presence  of  turning  points,  cannot  be 
parameterized  in  the  form  (X(J),  J).  The  upper  part  of  the  manifold  denotes  the  stable 
solutions  and  the  lower  part  the  unstable  ones,  assuming  there  are  no  Hopf  bifurcations. 

To  generate  this  solution  set,  (X,  J)  is  reparameterized  into  (Af(s),  J(s))  where  s  is  a 
new  independent  parameter  and  J  becomes  an  eigenvalue.  The  system  in  (11)  can  now  be 
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written 


J^(X(s),J(s))=a,  (12) 

and  the  dependence  of  s  on  the  augmented  solution  vector  {X,  J)  is  specified  by  an  extra 
scalar  equation 

Ar(X(s),J(s),s)  =  0,  (13) 

which  is  chosen  such  that  s  approximates  the  arclength  of  the  solution  branch  in  a  given 
phase  space  [16-17].  Rather  than  solving  the  coupled  system  in  (12)  and  (13),  we  replace  J 
by  a  function  J  oi  z  and  we  let  the  unknown  Z  =  (A*^,  JY  be  the  solution  of  a  three-point 
limit  value  problem 

H(Z,s)=  ^  =0,  (14) 

dz 

VM{X{z),J{zIs)\ 

where  z  is  a  given  point  in  [0,  L].  The  main  advantage  of  considering  (14)  is  that,  if  (12)  can 
be  discretized  in  a  block-tridiagonal  form,  then  (14)  has  the  same  property.  The  numerical 
Jacobian  matrix  obtained  from  (14)  has  a  block-tridiagonal  structure  and  one  can  prove 
easily  that  it  is  nonsingular  at  simple  turning  points.  Specially  developed  block  tridiagonal 
linear  equation  solvers  can  then  be  used  to  invert  the  corresponding  linear  systems. 

The  system  in  (14)  is  solved  by  combining  a  first-order  Euler  predictor  and  a  corrector 
step  involving  Newton-Uke  iterations  and  adaptive  gridding.  To  resolve  the  high  activity 
regions  of  the  dependent  solution  components,  the  mesh  must  be  refined  in  the  continuation 
calculations.  Specifically,  after  solving  the  governing  equations  on  a  previously  determined 
mesh,  we  determine  a  new  equidistributed  mesh 

Z‘  =  {0  =  zl<z;<...<z:.=L),  (15) 

by  imposing  the  condition 

dz  =  l,  (16) 

Jzr 

where  the  weight  function  w  depends  upon  the  gradient  and  curvature  of  the  dependent 
solution  components  and  on  mesh  regularity  properties.  The  governing  equations  are  solved 
on  this  new  mesh  and  then  another  continuation  step  is  taken.  It  optimizes  the  number  of 
points  and  determines  the  new  grid  in  only  one  pass.  In  addition,  the  method  also  takes 
into  account  every  component  of  the  solution  in  forming  a  new  mesh. 

Results  and  Discussion 

In  this  section  we  discuss  the  results  of  a  sequence  of  numerical  calculations  in  which 
the  dependence  of  the  soot  volume  fraction  as  a  function  of  the  strain  rate  is  investigated 
for  methane-air  and  ethylene-air  diffusion  flames.  All  the  calculations  were  done  at  one 
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atmosphere  with  the  fuel  and  oxidizer  streams  having  temperatures  of  298K.  The  lengths 
of  the  computational  domains  were  5.00  cm  for  ethylene  and  4.25  cm  for  methane.  A 
total  of  20  soot  size  classes  was  used  in  the  computations.  Each  continuation  calculation 
contained  between  125-175  adaptively  determined  grid  points  and  all  of  the  calculations 
were  performed  on  an  IBM  RS/6000  Model  560  computer.  Approximately  1000  solutions 
taking  75  hours  of  CPU  time  were  needed  to  generate  the  upper  branch  of  the  C-shape 
solution  curve  for  ethylene  and  approximately  350  solutions  taking  25  hours  of  CPU  time 
were  needed  for  methane. 

In  Figure  4  we  display  the  variation  of  the  peak  temperature  vs.  the  strain  rate  for 
both  sets  of  flames.  Notice  the  large  difference  in  the  extinction  strain  rates  between  the 
two  cases.  Using  the  strain  rate  defined  in  (10),  we  determined  an  extinction  strain  rate 
of  338  sec-^  for  methane  and  1150  sec'^  for  ethylene.  When  the  strain  rate  is  defined  as 
the  maximum  axial  velocity  gradient  on  the  oxidizer  side,  we  compute  an  extinction  strain 
rate  for  methane  of  620  sec'^  With  this  same  strain  rate  definition,  the  experimental 
value  reported  in  [6]  is  380  sec'^  Since  this  computed  value  is  far  outside  what  should 
be  considered  acceptable  agreement,  the  strain  rate  at  extinction  has  also  been  computed 
using  a  mechanism  utilized  in  our  prior  study  (Hall,  Smooke,  and  Colket,  1997).  With  that 
mechanism,  the  strain  rate  at  extinction  was  calculated  to  be  390  sec  ^  The  peak  volume 
fractions  as  predicted  from  this  earUer  mechanism  at  low  strain  rates  were  comparable  to 
those  from  the  modified  version  of  GRIMech. 

Calculations  neglecting  radiation  losses  were  also  performed.  Radiation  effects  are 
relatively  small  for  ethylene  flames  for  strain  rates  above  approximately  100  sec'^  (about 
13K  depression  at  40  sec-^  and  6K  depression  at  80  sec"^).  For  lower  strain  rates  the  gas 
cooUng  effects  of  radiation  are  important  and  have  a  significant  influence  on  both  methane 
and  ethylene  temperature  profiles.  In  particular,  we  have  computed  methane  flames  down 
to  strain  rates  of  6  sec"^;  for  these  low  strain  rate  flames  radiative  cooling  effects  ^e  on 
the  order  of  lOOK.  For  ethylene  flames,  calculations  were  calculated  down  to  15  sec"^;  the 
corresponding  cooling  effect  is  45K.  The  rapid  reduction  in  the  peah  temperature  for  the 
both  the  ethylene  and  methane  flames  at  low  strain  may  be  unphysical  due  to  the  neglect 
of  self-absorption  inherent  in  the  optically  thin  model. 

The  corresponding  profiles  of  maximum  soot  volume  fraction  are  shown  in  Figure  5  for 
the  methane-air  and  ethylene-air  flames.  The  very  low  values  for  the  methane  flame  and  its 
rapid  decline  are  consistent  with  the  known  difficulty  of  stabilizing  a  sooting  methane-air 
flame  in  the  laboratory.  The  variation  of  volume  fraction  in  the  ethylene  flames  has  b^n 
shown  to  be  reasonably  consistent  with  experimental  measurements  of  volume  fraction 
strain  rate  sensitivity  [25]. 

The  importance  of  fully  coupling  the  radiation  and  scrubbing  processes  within  the 
flame  and  soot  equations  can  be  clearly  demonstrated  under  low  strain  rate  conditions. 
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The  45K  increase  in  temperature  when  thermal  radiation  loss  was  excluded  in  ethylene  at 
15  sec"^  increases  peak  volume  fraction  by  75%.  Furthermore,  when  both  scrubbing  and 
radiation  effects  are  not  included  in  the  model  the  temperature  increases  by  another  25K 
and  soot  by  another  factor  of  two.  These  effects  are  even  more  dramatic  in  the  case  of  the 
low  strain  rate  methane  flame  at  6  sec"^  where  the  lOOK  increase  in  flame  temperature 
with  radiation  excluded  leads  to  a  factor  of  five  increase  in  the  peak  soot  volume  fraction! 
As  stated,  however,  the  neglect  of  self-absorption  undoubtedly  overestimates  the  radiation 
losses  and,  thus,  the  predictions  in  Fig.  5  are  likely  low. 

Figures  6  and  7  show  the  variation  in  peak  mole  fractions  of  key  growth  species  for  the 
ethylene  and  methane  flames.  Of  interest  are  the  relatively  ’flat’  profiles  of  the  acetylene 
concentration  (a  fantor  of  three  change)  and  of  the  benzene  concentration  (only  an  order 
of  magnitude  decrease).  These  numbers  are  contrasted  with  the  very  dramatic  drop  in  the 
soot  volume  fraction  nearly  five  orders  of  magnitude  decay  for  the  ethylene  flame  (Figure 
5).  While  hydrogen  is  also  decaying,  the  change  in  species  concentrations  alone  cannot 
account  for  the  dramatic  drop  in  soot  volume  fraction.  The  temperature  dependence  of  the 
inception  process  is  small  (10-13  kcal/mole),  while  that  of  the  growth  rate  is  noticeably 
greater  (31.8  kcal/mole).  The  change  in  this  latter  rate  is  only  about  a  factor  of  three  for 
the  calculated  temperature  range  (1750-2100K).  While  each  of  these  processes  contributes 
to  the  large  decay  in  soot  volume  fraction,  together  they  cannot  account  for  the  five  order 
of  magnitude  decay.  As  shown  in  [11],  oxidation  plays  a  relatively  minor  role  in  soot 
production  in  opposed  jet  flames;  hence,  an  increase  in  the  oxidation  rate  probably  does 
not  explain  the  decline  either.  A  remaining  explanation  is  the  very  significant  decrease  in 
residence  time.  For  the  strain  rate  variation  of  15  sec“^  to  1150  sec“‘  for  the  ethylene  flame, 
this  represents  approximately  a  factor  of  75  decrease  in  residence  time.  This  parameter 
is  likely  the  most  important  factor  in  the  significant  decay  of  soot  volume  fraction  in  the 
ethylene  flame;  and  the  decay  in  the  concentration  of  the  first  aromatic  ring  is  the  second 
most  important  parameter. 

Acetylene  levels  in  the  methane  flame  (Figure  7)  are  slightly  lower  than  in  the  ethylene 
flame  and  benzene  is  about  an  order  of  magnitude  lower  (Figure  6).  These  differences 
as  well  as  the  lower  peak  flame  temperatures  contribute  to  the  much  lower  soot  from  the 
methane  flame. 

Figure  8  shows  a  comparison  of  the  soot  variation  with  strain  rate  for  the  present 
model  and  that  of  the  model  of  Fairweather  et  al.  [29].  The  inception  rate  and  the 
surface  growth  rate  in  the  latter  model  are  both  proportional  to  acetylene  with  activation 
energies  of  42.0  and  24.0  kcal/mole,  respectively.  This  is  usually  a  monodisperse  model; 
in  this  work,  it  has  been  incorporated  into  the  sectional  representation.  The  faet  that  the 
relative  variation  with  strain  rate  is  similar  for  the  two  models,  in  spite  of  the  considerable 
functional  differences  and  the  fact  that  the  Fairweather,  et  al.  model  is  dependent  only 
on  the  acetylene  concentration,  strongly  supports  the  previous  agrument  regarding  the 
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importance  of  residence  time.  The  fact  that  peah  soot  volume  fractions  predicted  from 
our  model  fall  off  slightly  more  rapidly  than  those  predicted  using  the  Fairweather  model 
is  presumably  due  to  a  combination  of  the  lower  activation  energy  of  the  surface  growth 
step  in  the  Fairweather  model  and  the  order  of  magnitude  decrease  in  benzene  to  which 
our  inception  model  is  sensitive. 

Conclusions 

Arclength  continuation  methods  have  been  incorporated  into  an  opposed  jet  soot  pre¬ 
diction  code  with  complex  chemistry,  detailed  particle  transport/kinetics,  and  radiation. 
The  continuation  technique  greatly  reduces  the  effort  required  to  carry  out  parametric 
variations  with  opposed  jet  programs.  Calculations  have  been  earned  out  for  the  strain 
rate  dependence  of  atmospheric  pressure  methane-air  and  ethylene-air  coimterflow  diffu¬ 
sion  flames.  Comparisons  with  experimental  data,  examination  of  the  effects  of  thermal 
radiation,  and  discussion  of  the  main  factors  governing  the  strain  rate  dependence  of  soot 
formation  in  these  flames  have  been  presented. 
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Figure  Captions 

1.  Coordinate  system  for  counterflow  flame  analysis. 

2.  Comparison  of  computed  and  experimental  velocity  profiles.  Data  for  opposed  jet  flames 
#1  and  #4  of  Wang  et  ad.  [25]. 

3.  Comparison  of  computed  and  experimental  soot  volume  fraction  profiles.  Data  for 
opposed  jet  flames  ^1  and  ^4  of  Wang  et  ad.  [25]. 

4.  Peak  temperature  variation  as  a  function  of  strain  rate  for  ethylene  and  methane-air 
counterflow  diffusion  flames. 

5.  Peak  volume  fraction  variation  as  a  function  of  strain  rate  for  ethylene  and  methane-air 
counterflow  diffusion  flames. 

6.  Variation  of  acetylene  and  benzene  mole  fractions  with  strain  rate  for  ethylene-air 
counterflow  diffusion  flames. 

7.  Variation  of  acetylene  and  benzene  mole  fractions  with  strain  rate  for  methane-air 
counterflow  diffusion  flames. 

8.  Comparison  of  peak  soot  volume  fraction  versus  strain  rate  for  present  model  and  the 
model  of  Fairweather  et  al.  [29]. 
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Abstract 

A  detailed  soot  growth  model  in  which  the  dynamical  equations  for  particle  produc¬ 
tion  have  been  coupled  to  the  flow  and  gaseous  species  conservation  equations  has  been 
developed  for  an  axisymmetric,  laminar,  coflow  diffusion  flame.  Results  from  the  model 
have  been  compared  to  experimental  data  for  a  confined  methane-air  flame.  The  two- 
dimensional  system  couples  detailed  transport  and  flnite  rate  chemistry  in  the  gas  phase 
with  the  aerosol  dynamical  equations  in  the  sectional  representation.  The  formulation 
includes  detailed  treatment  of  the  transport,  inception,  surface  growth,  oxidation,  and 
coalescencfe  of  soot  particulates.  Effects  of  thermal  radiation  and  particle  scrubbing  of 
gas  phase  growth  and  oxidation  species  are  also  included.  Predictions  and  measurements 
of  temperature,  soot  volume  fractions  and  selected  species  are  compared  over  a  range  of 
heights  and  as  a  function  of  radius.  Flame  heights  are  slightly  overpredicted  and  local 
temperatures  and  volume  fractions  are  underpredicted.  We  believe  the  inability  to  repro¬ 
duce  accurately  bulk  flame  parameters  directly  inhibits  the  ability  to  predict  soot  volume 
fractions  and  these  differences  are  likely  a  result  of  uncertainties  in  the  experimental  inlet 
conditions.  Predictions  of  the  distributions  of  particle  sizes  indicate  the  existence  of  (rela¬ 
tively)  low  molecTilar  weight  species  along  the  centerline  of  the  burner  and  trace  amoimts 
of  the  particles  which  escape  from  the  flame,  unoxidized.  Oxidation  of  particulates  is  dom¬ 
inated  by  reactions  with  hydroxyl  radicals  which  attain  levels  approximately  ten  times 
higher  than  calculated  equilibrium  levels.  Gas  cooling  effects  due  to  radiative  loss  are 
shown  to  have  a  very  significant  effect  on  predicted  soot  concentrations. 
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COMPUTATIONAL  AND  EXPERIMENTAL  STUDY  OF 
SOOT  FORMATION  IN  A  COFLOW,  LAMINAR  DIFFUSION  FLAME 


1.  Introduction 

In  the  last  few  years,  we  have  witnessed  a  shift  in  combustion  research.  The  push  for 
higher  combustion  efficiency  in  propulsion  applications  that  dominated  much  of  research 
in  the  past  few  decades  is  gradually  being  replaced  by  a  drive  towards  cleaner  combustion. 
This  is  a  direct  result  of  environmental  consciousness  and  it  has  been  translated  into  stricter 
air  quality  legislation.  Although  it  originated  as  a  reaction  to  regulatory  pressure,  research 
in  the  area  of  pollutant  formation  and  control  will  become  economically  indispensable  to 
the  export  of  combustion-related  technologies  and  products  worldwide.  In  particular,  as 
emissions  legislation  becomes  more  restrictive,  a  detailed  understanding  of  soot  formation 
and  unburned  hydrocarbons  in  flames  will  be  critical  for  the  design  of  pollutant  abatement 
strategies  and  for  the  preservation  of  the  competitiveness  of  combustion  related  industries. 

Modeling  soot  formation  in  practical  combustion  systems  is  an  extremely  challenging 
problem.  Reacting  flow  problems  for  simple  fuels  such  as  methane  can  require  the  solution 
of  more  than  50  chemical  species  in  addition  to  the  temperature  and  the  fluid  dynamic 
variables.  The  inclusion  of  soot  inception,  growth  and  oxidation  processes  in  these  models 
greatly  increases  the  level  of  complexity  of  the  problem.  In  fact,  due  to  the  number  of 
dependent  ,imknowns  and  the  resolution  with  which  reaction  fronts  need  to  be  resolved, 
gas  phase  Combustion  calculations  with  detailed  chemistry  have  only  recently  moved  from 
one-  to  two-  and  some  three-dimensional  systems.  The  increased  effort  with  which  mul¬ 
tidimensional  combustion  phenomena  have  been  studied  numerically  is  attributable  in 
part  to  advances  in  the  development  of  computational  algorithms,  and  in  part  to  higher 
speed/larger  memory  workstations.  This  combination  has  enabled  the  combustion  scien¬ 
tist  to  investigate  chemically  reacting  systems  that  were  computationally  infeasible  only  a 
few  years  ago.  Nevertheless,  the  investigation  of  soot  formation  with  detailed  chemistry 
in  a  generic  multidimensional  configuration  is  still  beyond  our  current  computational  abil¬ 
ity.  The  laminar  diffusion  flame,  however,  provides  a  natural  environment  in  which  one 
can  investigate  the  interaction  of  soot  formation  with  detailed  gas-phase  chemistry  in  a 
multidimensional  system. 

Most  studies  using  detailed  chemical  kinetics  and  coupled  models  of  soot  production 
and  oxidation  have  focused  on  one-dimensional  geometries.  A  few  recent  studies  [Kaplan 
et  aJ.,  1996,  Kennedy  et  al.,  1996]  have  focused  on  jet  diffusion  flames,  using  simplified, 
monodisperse,  soot  formation  models  with  skeletal  kinetic  mechanisms.  In  this  paper  we 
modify  the  sectional  soot  formation  model  developed  in  Colket  and  Hall  (1994)  and  Hall  et 
al.  (1997)  for  incorporation  into  a  laminar,  axisymmetric,  diffusion  flame  in  which  a  cylin- 
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drical  fuel  stream  is  surrounded  by  a  coflowing  oxidizer  jet.  Computationally,  we  employ 
a  velocity- vorticity  model  (Ern  et  al.,  1995)  in  which  the  governing  conservation  equations 
are  solved  with  detailed  transport  and  finite  rate  chemistry  submodels  to  predict  the  tem- 
peratme,  species  mass  fractions  and  velocity  fields  as  functions  of  the  two  independent 
coordinates.  Appropriate  sectional  equations  are  included  for  describing  the  different  soot 
particle  size  classes  with  a  convective,  diffusive  and  chemical  production  balance.  These 
equations,  along  with  surface  growth,  oxidation  and  radiation  (from  both  gas  and  par¬ 
ticulate  species)  are  fully  integrated  with  the  governing  equations.  A  discrete  solution  is 
obtained  on  a  two-dimensional  grid  by  employing  Newton’s  method  with  adaptive  mesh 
refinement.  In  many  models  diffusion  in  the  axial  direction  is  neglected:  here,  we  consider 
the  fully  elliptic  problem.  Experimentally,  we  utilize  thermocouple  particle  densitometry 
and  mass  spectrometry  to  obtain  profiles  of  the  soot  volume  fraction  and  the  species  con¬ 
centrations  at  various  heights  in  the  flame  (McEnally  and  Pfefferle,  1997).  We  apply  the 
computational  model  and  experimental  methods  to  a  lightly  sooting  methane-air  diffusion 
flame.  In  the  next  section,  we  formulate  the  problem.  In  Section  3,  we  describe  the  exper¬ 
imental  approach.  The  computational  method  is  discussed  in  Section  4  and  results  of  the 
study  are  contained  in  Section  5. 

2.  Problem  Formulation 

While  the  ultimate  goal  in  combustion  modeling  is  the  solution  of  three-dimensional 
turbulent  reacting  flows  with  finite  rate  chemistry,  there  are  still  important  less  complex 
systems  that  can  be  analyzed  in  detail  with  current  computational  resources.  The  ax- 
isymmetric,  laminar,  diffusion  flame  is  one  such  configuration  (Smooke  et  al.,  1990,  1992, 
1996;  Xu  et  al.,  1993;  Ern  et  al.,  1995)  This  is  the  flame  type  of  most  practical  combus¬ 
tion  devices  and  it  is  important  in  understanding  interactions  of  heat  and  mass  transfer 
with  chemical  reactions  in  gas  turbines  and  commercial  burners.  The  ability  to  predict 
the  coupled  effects  of  complex  transport  phenomena  with  detailed  chemical  kinetics  in 
these  systems  is  critical  in  the  modeling  of  turbulent  reacting  flows,  in  improving  engine 
efficiency  and  in  understanding  the  processes  by  which  pollutants  are  formed. 

2.1  Soot  Modeling 

Soot  kinetics  axe  modeled  as  coalescing,  solid  carbon  spheroids  undergoing  surface 
growth  in  the  free  molecule  limit.  The  particle  mass  range  of  interest  is  divided  into 
sections  (Gelbard  and  Seinfeld,  1980),  and  a  dynamical  equation  is  written  for  each  section 
including  coalescence,  surface  growth,  and  oxidation.  Sectional  analysis  makes  it  possible 
to  obtmn  the  particle  size  distribution  without  a-priori  assumptions  about  the  form  of  the 
distribution.  For  the  first  bin,  an  inception  source  term  is  included.  The  incorporation  of 
these  dynamical  equations  into  a  transport /conservation  equation  for  each  section  includes 
both  thermophoresis  and  an  effective  bin  diffusion  rate.  In  the  gas  species  conservation 
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equations,  provision  is  made  for  complex  chemistry;  additional  source  terms  are  included 
to  account  for  scrubbing  or  generation  of  gaseous  species  arising  from  the  particle  growth 
and  oxidation  processes.  The  gas  and  soot  are  additionally  coupled  through  non-adiabatic 
radiative  loss  from  both  the  gas  and  the  soot  in  the  optically-thin  approximation.  For  one 
atmosphere  calculations,  the  sectional  coefficients  for  transport,  coalescence  and  surface 
processes  are  calculated  in  the  free  molecule  limit  (Gelbard  and  Seinfeld,  1980),  with 
surface  growth  and  oxidation  proportional  to  particle  surface  area.  In  the  free  molecule 
limit,  the  particle  mass  and  medium  dependencies  of  the  sectional  coefficients  are  factored 
in  such  a  way  that  they  only  have  to  be  evaluated  once.  Oxidation  of  soot  by  O2  and 
OH  is  treated  as  described  in  (Hall  et  al.,  1997).  The  inception  model  employed  here  is 
based  on  an  estimate  of  the  formation  rate  of  two-  and  three-ringed  aromatic  species,  and 
is  a  function  of  local  acetylene,  benzene,  phenyl  and  molecular  hydrogen  concentrations. 
Assuming  steady-state  values  of  intermediates,  and  also  that  H2  >>  C2H2,  the  rate  of 
production  of  the  polyaromatic  species  can  be  estimated  to  be  (Hall  et  ai.,  1997) 

=  10"  cc/mole/sec,  (2.1) 


and 


cc/mole/sec,  (2.2) 

/ 

where  the  gas  phase  concentrations  and  temperatures  are  evaluated  at  local  conditions. 
With  the  further  assumptions  that  inception  is  limited  by  the  formation  of  polyaromatics 
and  oxidation/decomposition  of  such  species  can  be  neglected,  the  inception  rate,  Si,  in 
grams/cc/sec,  was  initially  assigned  to 


S,  =  127  X  3%^  +  178  X  (2.3) 

at  at 

where  the  constants  (molecular  weights)  are  provided  to  convert  from  molar  to  mass  units. 
The  second  term  in  equation  2.3  plays  a  negligible  role  in  methane-air  diffusion  flames  but 
has  been  included  here  for  generality.  The  contributions  from  both  inception  processes  are 
incorporated  in  the  first  sectional  bin,  whose  lower  mass  boundary  is  set  equal  to  the  mass 
of  the  smallest  inception  species  (See  Section  5).  In  the  sectional  representation  (Gelbard 
and  Seinfeld,  1980),  the  sectional  mass  boundaries  vary  linearly  on  a  logarithmic  scale. 

The  surface  growth  rate  is  nominally  that  of  Harris  and  Weiner  (1983),  who  report  a 
rate  proportional  to  acetylene  concentration;  their  rate  is  normalized  in  this  application 
to  1650K,  with  a  31.8  kcal  activation  energy  as  suggested  by  Hura  and  Glassman  (1988). 
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It  has  been  found  (Smooke  et  aJ.,  1997)  that  enhancing  the  Harris  and  Weiner  rate  by  a 
factor  of  two  results  in  predictions  of  soot  levels  in  opposed  jet  ethylene  flames  in  good 
agreement  with  measurements  of  Wang  et  aJ.  (1997).  The  enhancement  by  a  factor  of  two 
brings  the  Harris  and  Weiner  rate,  determined  from  studies  on  premixed  ethylene  flames, 
into  better  agreement  with  results  obtained  in  diffusion  flames  (Axelbaum  et  aJ.,  1988, 
Kennedy  et  ad.,  1990,  Sunderland  and  Faeth,  1996). 

The  surface  growth/inception  models  used  in  this  paper  agree  well  with  available  ethy¬ 
lene  opposed  jet  data.  Nevertheless,  it  is  not  suggested  that  they  represent  a  universal 
solution  to  the  problem  of  soot  formation  because  comparisons  with  soot  data  in  other 
fuels  have  not  been  carried  out.  The  possibility  of  aging  of  the  surface  growth  rate  has 
not  been  accounted  for  here,  nor  are  the  results  with  finite  coalescence  appropriate  for 
aggregate  formation.  These  effects  should  be  included  in  future  modeling  efforts.  Fun¬ 
damental  questions  relating  to  inception,  temperature  dependence  of  surface  growth  and 
precise  dependence  of  surface  growth  on  gas  phase  species  concentrations  remain.  Whether 
soot  growth  rates  derived  from  premixed  and  opposed  jet  studies  are  self-consistent  is  still 
an  open  issue  as  weU. 

A  primary  parameter  of  interest  in  this  paper  is  the  soot  volume  fraction,  for  which 
convergence  is  generally  obtained  with  a  smaller  number  of  sections  than  is  the  case  for 
particle  size  and  number  density  (Hall  et  ad.,  1997).  The  number  of  sections  required  for 
convergence  must  be  examined  in  each  case,  although  evidence  points  to  the  relative  magni¬ 
tudes  of  surface  growth  and  inception  as  the  important  parameter.  When  the  contribution 
of  inceptiop  to  total  volume  fraction  is  significant,  although  not  necessarily  dominating, 
the  numbdr  of  sections  required  can  be  very  small  (Hall  et  ad.,  1997),  with  larger  numbers 
necessary  when  the  contribution  of  surface  growth  to  volume  fraction  becomes  very  much 
larger  than  that  from  inception. 

2.2  Governing  Equations 

If  we  combine  the  gas-phase  diffusion  flame  model  which  employs  a  velocity- vorticity 
formulation  (Ern  et  ad.,  1995)  with  the  sectional  approach  presented  in  Hall  et  ad.  (1997), 
the  coupled  set  of  governing  partial  differential  equations  can  be  written  in  the  form 

Radial  Velocity 

d^Vr  d^Vr 
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Axial  Velocity 


d^v,  doj  d  ^v^dp\  d  /Vr\  d  fvrdp\ 

dr^  dz^  dz  dz  \  p  dz  j  dz  \  r  )  dz  \  p  dr )  ' 

Vorticity  Transport 


d  fu)  dtp  \  d  P  CO  dtp  '\ 
dz  \r  dr  j  dr  \r  dz  j 


=  r 


dp 

-^^Tr 


rV 


•  iso(/>) 
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(2.5) 


(2.6) 


Species 


pv, 


’■  dr 


Wtiwk+wl)  =  (i,  k  =  l,2,...,K, 


(2.7) 


Sections 


dY,  dY,  Id 
+  r-Fr 


(2.8) 


Energy 


dz 


K 


+  +  Wl)  —  V  •  =  0, 

ib=l 


(2.9) 


Equation  of  State 


pW 

RT‘ 


(2.10) 
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The  system  is  closed  with  appropriate  boundary  conditions  on  each  side  of  the  computa¬ 
tional  domain.  For  a  confined  coflowing  flame  we  have 

Axis  of  Symmetry  (r  =  0): 


dvg  _  du} 
dr  dr 

Exit  (z  — »  oo): 


^  =  ^  =  n 

dr  dr  ’ 


k  =  +  M, 


dVi  _  doJr 

IT 


dY, 

dz 


dz 


=  0, 


k  =  1,2,...,  K  -|- 


Inlet  {z  =  0): 
r  <  Ri 


Vr  =  0, 


Rf  <r  <  Ro 


/ 


Vz  =  Vf, 
dVr  dVt 

n  =  k  =  l,2,..^,K  +  M, 
T  =  Tf, 

Vr  =  0, 


(2.11) 


(2.12) 


(2.13) 


Wall  (r  =.iio): 


Vz  =  VoXi 
dvr  dvz 

n  =  ^  =  i,2,...,/i:+M, 

T  =  Tox. 


(2.14) 


dr 


Vr=Vz=  0, 

dvz 

0,  k  =  1,2, . . .  1 K  M, 
T  =  T^uatl- 


(2.15) 


The  subscripts  I  and  O  refer  to  the  inner  jet  and  the  outer  jet,  respectively,  and  Vf,voxi 
Ykp  1  Ykox  1  ^F,  Tox  and  T,ji,aii  are  specified  quantities. 
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In  addition  to  the  variables  already  defined,  K  represents  the  total  number  of  gas 
phase  species  and  M,  the  number  of  soot  (size)  sections;  T,  denotes  the  temperature;  F*,, 
the  mass  fraction  of  the  fc"*  gas-phase  species  {k  <  K)  or  the  (k  —  K)*''  soot  size  class 
(k  >  K)\  p,  the  pressure;  Vr  and  v^,  the  velocities  of  the  fluid  mixture  in  the  radial  and 
axial  directions,  respectively;  p,  the  mass  density;  Wt,  the  molecular  weight  of  the  fc*'* 
species;  W,  the  mean  molecular  weight  of  the  mixture;  J2,  the  universal  gas  constant;  A, 
the  thermal  conductivity  of  the  mixture;  Cp,  the  constant  pressure  heat  capacity  of  the 
mixture;  Cpk,  the  constant  pressure  heat  capacity  of  the  k*''  species;  tbt,  the  molar  rate 
of  production  of  the  k*^  species  per  unit  volume;  w^,  the  molar  rate  of  production  of  the 
k*'*  species  per  unit  volume  due  to  scrubbing/replenishment  by  the  soot  growth/oxidation 
processes;  the  rate  of  change  of  section  k  due  to  inception,  surface  growth,  oxidation 
and  coalescence;  hk,  the  specific  enthalpy  of  the  k*^  species;  g,  the  gravitational  constant; 
p,  the  viscosity  of  the  mixture;  and  14,,  the  diffusion  velocities  of  the  k^’'  species  in 
the  radial  and  axial  directions,  respectively;  V^r,  Vtz,  14%  14%  sectional  thermophoretic  and 
diffusion  velocities  of  the  k*^  soot  size  class  and  V  •  qr  the  divergence  of  the  net  radiative 
flux  for  gas  bands  and  soot  in  the  optically-thin  limit. 

We  write  the  diffusion  velocities  in  the  r  and  z  directions  in  the  form 


(2.16) 

Vk^  =  -(l/Xk)Dk^,  k  =  l,2,.. 

(2.17) 

where  Xk  is  the  mole  fraction  of  the  fc*'*  species  and  Dk  is  related  to  the  binary  diffusion 
coefficients'' through  the  relation  (see,  e.g.,  Curtiss  and  Hirschfelder,  1949) 

(1  -  n) 

The  sectional  thermophoretic  velocities  in  the  free  molecule  regime  are  given  by  (Hall  et 
al.,  1997) 

f 

The  sectional  diffusion  velocities  are  written  as  in  (2.16-2.17)  with  a  mass  weighted  mean 
diffusion  coefficient  for  each  bin  (Hall  et  aJ.,  1997).  The  binary  diffusion  coefficients,  the 
viscosity,  the  thermal  conductivity  of  the  mixture,  the  chemical  production  rates  and  the 
thermodynamic  quantities  were  evaluated  using  vectorized  and  highly  optimized  transport 
and  chemistry  libraries  (Giovangigli  and  Darabiha,  1987). 

Anticipating  that  radiative  losses  could  have  a  significant  influence  on  soot  levels  com¬ 
pared  to  a  nonradiative  system,  we  have  included  an  optically-thin  radiation  model  in 
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our  calculations  (Hall,  1993,  1994).  We  assume  that  for  methane-air  mixtures  the  only 
significant  radiating  species  are  H^O,  CO  and  CO 2-  By  utilizing  an  optically- thin  limit  in 
which  self-absorption  of  radiation  is  neglected,  the  divergence  of  the  net  radiative  flux  can 


be  written 

V-qR  =  Cf„T^  +  47r  otikPkIi.,, ,  (2.21) 

ik 


where  /„  is  the  soot  volume  fraction  and  is  the  Planck  function  evaluated  at  the 
gas  band  centers  of  the  contributing  vibration-rotation  or  pure  rotational  bands  whose 
integrated  intensities  are  given  by  Oij,. 

3.  Experimental  Methods 

3.1  Burner  and  Reactants 

Atmospheric  pressure,  over  ventilated,  axisymmetric,  coflowing,  nonpremixed  laminar 
flames  were  generated  with  a  burner  in  which  the  fuel  flows  from  an  uncooled  12  mm  diam¬ 
eter  vertical  brass  tube  and  the  oxidizer  flows  from  the  annular  region  between  this  tube 
and  a  102  inm  diameter  concentric  tube  (see  Figure  1).  Details  of  the  burner  construc¬ 
tion  and  operation  are  given  elsewhere  (McEnally  and  Pfefferle,  1997).  The  oxidizer  was 
air  (710  cm®/min),  while  the  fuel  was  a  mixture  containing  methane  (240  cm^/min)  and 
argon  (2.4  cm®/min).  The  added  argon  edded  species  measurements  as  discussed  below. 
All  flowrates  were  governed  by  electronic  mass  flow  controllers  and  are  estimated  to  be 

accurate  to  within  5%. 

/ 

3.2  Temperature  and  Soot  Volume  Fraction  Measurements 

Temperatures  were  measured  with  uncoated  75  pm.  wire-diameter  type  R  thermocou¬ 
ples.  By  measuring  profiles  in  coflowing  flames  with  various  size  thermocouples,  we  have 
determined  that  while  conduction  errors  of  up  to  lOOK  occur  at  some  flame  locations  with 
125/xm  wires,  they  are  negligible  with  75  pm  wires  (McEnally  et  aJ.,  1997).  The  measured 
junction  temperatures  were  converted  to  gas  temperatures  with  a  heat  balance  (radiation 
energy  loss  =  convection  energy  gain)  and  the  following  parameters:  Nusselt  number  = 
2.26,  junction  diameter  =  240  pm  (measured  with  a  microscope),  temperature-dependent 
thermal  conductivity  data  for  air  (Holman,  1986),  and  temperature-dependent  type  S 
thermocouple  emissivities  (Bradley,  1961).  We  estimate  that  the  absolute  error  in  the  gas 
temperatures  is  less  than  50K  in  soot-free  flame  regions.  In  soot-containing  flame  regions, 
particles  rapidly  deposit  onto  the  junction,  rendering  the  junction  diameter  and  emissivity 
values  listed  above  incorrect.  Therefore  the  junction  was  cleaned  before  every  measure¬ 
ment  by  retracting  it  into  the  oxidizing  layer  outside  the  flame,  and  the  temperature  was 
recorded  quickly  (0.75  seconds)  after  the  junction  was  reinserted  into  the  flame. 

Soot  volume  fractions  were  measured  with  a  new  technique,  thermocouple  particle 
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densitometry  (TPD),  that  was  suggested  by  Eisner  and  Rosner  (1985)  and  recently  imple¬ 
mented  by  McEnally  et  aJ.  (1997).  A  clean  thermocouple  was  introduced  into  the  flame 
and  its  junction  temperature  was  recorded  periodically  for  200  seconds;  the  soot  volume 
fraction  was  then  determined  by  optimizing  the  agreement  between  the  measured  history 
and  one  calculated  with  the  thermophoretic  mass  transfer  formulation  of  Eisner  and  Ros¬ 
ner  (1985).  In  the  work  cited  above  we  demonstrated  that  TPD  measurements  agreed 
with  soot  volume  fractions  measured  by  the  conventionzil  technique  of  laser  extinction  in 
counterflow  and  coflowing  nonpremixed  flames.  We  developed  TPD  and  employed  it  here 
because  the  combination  of  low  soot  volume  fractions  (<1  ppm)  and  short  pathlengths  (<5 
mm)  make  laser  extinction  measurements  difficult  in  methane  coflowing  flames  (Garo  et 
al.,  1986),  (Shaddix  et  aJ.,  1994).  With  TPD,  soot  volume  fractions  can  be  easily  measured 
down  to  at  least  0.02  ppm  by  allowing  sufficient  time  to  pass  for  a  measurable  change  in  the 
junction  temperature  to  occur.  The  technique  also  allows  spatially-resolved  measurements 
without  tomographic  reconstruction.  The  results  have  an  absolute  accuracy  of  ±  50%  and 
a  relative  accuracy  of  ±  10%. 

Recent  experiments  have  shown  that  nonpremixed  ffames  contain  two  kinds  of  soot 
particles:  an  early  type  that  appears  translucent  when  viewed  under  an  electron  micro¬ 
scope,  and  a  later  type  that  appears  dark  (Mergaridis  and  Dobbins,  1989).  Similarly,  two 
types  of  particles  have  been  identified  in  premixed  flames,  where  the  early  particles  were 
shown  to  be  transparent  to  visible  light  (d’Alessio  et  al.,  1992).  Results  from  our  labo¬ 
ratory  suggest  that  soot  volume  fractions  measured  by  TPD  include  both  particle  types, 
while  laser  extinction  and  laser-induced  incandescence  respond  to  only  the  dark  particles 
(McEnally /and  Pfefferle,  1996),  (Koylii  et  al.,  1996). 

3.3  Species  Measurements 

Species  concentrations  were  measured  by  extracting  gases  from  the  flames  with  a 
narrow-tip  quartz  microprobe,  and  analyzing  them  with  on-line  mass  spectrometry.  The 
sampling  procedure  and  its  performance  are  discussed  in  greater  detail  in  McEnally  and 
Pfefferle  (1997)  where  it  was  concluded  that  the  microprobe  spatial  resolution  was  roughly 
one  millimeter,  reactions  in  the  microprobe  tip  have  only  a  minor  affect  on  the  measured 
concentrations,  and  condensation/adsorption  losses  to  the  sample  manifold  walls  were  in¬ 
significant  for  biphenyl  (C12H10;  154  AMU)  and  all  lighter  species.  Following  Puri  (1992), 
we  have  equipped  the  probe  with  an  oscillating  sapphire  fiber  that  extends  through  its 
orifice  and  grinds  down  any  soot  deposits  that  accumulate  there.  This  prevents  complete 
clogging  of  the  orifice,  enabling  us  to  make  measurements  in  particle-containing  flame 
regions. 

Most  species  were  quantified  with  a  custom-built  photqionization/time-of-flight  mass 
spectrometer  (PTMS)  that  can  simultaneously  measure  all  C3  and  larger  hydrocarbons 
(except  propane  and  butane)  with  part-per-million  detection  limits  (Bermudez  and  Pf- 
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efFerle,  1995).  This  instrument  generates  ions  via  single-photon  photoionization  with  a 
low-intensity  vacuum-ultraviolet  (118.2  nm)  laser  beam.  The  calibration  factor  relating 
ion  signals  to  flame-gas  concentrations  was  determined  for  benzene  by  flowing  benzene/ air 
mixtures  of  known-composition  past  a  needle  valve  in  the  sample  manifold.  The  same 
calibration  factor  was  used  for  all  other  species  for  several  reasons:  independent  calibra¬ 
tions  would  have  been  difficult  since  the  detected  species  include  room-temperature  gases, 
liquids,  and  sohds;  several  masses  contain  multiple  isomers  whose  ratios  could  not  be  de¬ 
termined;  and  pure  samples  of  several  detected  species  could  not  be  easily  obtained.  We 
estimate  that  the  absolute  concentrations  of  benzene  are  accurate  to  30%  while  those  of 
other  species  are  accurate  to  within  a  factor  of  three.  The  primary  quantity  of  interest  here, 
the  relative  concentration  of  a  given  species  with  and  without  additives,  is  independent  of 
the  choice  of  calibration  factors  and  accurate  to  30%. 

A  commercial  electron-impact/quadrupole  mass  spectrometer  (EQMS;  Ley  bold  Inficon  . 
Transpector  CIOOM)  was  used  to  measure  methane,  acetylene,  and  argon.  The  EQMS 
ionizes  target  gases  with  a  fixed-energy  70  eV  electron  gun,  which  produces  extensive  ion 
fragmentation.  Methane  was  measured  with  its  CHJ  fragment-ion  to  avoid  interferences 
from  the  0+  fragment-ions  of  water  and  carbon  dioxide,  argon  with  its  Ar++  ion  to  avoid 
interferences  from  CaHJ  ions,  and  acetylene  with  its  parent  ion.  The  methane  and  acetylene 
concentrations  were  directly  calibrated  and  are  accurate  to  within  30%.  However,  ethylene 
fragment  ions  may  account  for  up  to  30%  of  the  acetylene  signal;  ethylene  and  its  C2HJ 
fragment-ion  were  completely  obscured  by  nitrogen  at  all  flame  locations,  so  only  an  upper 

botmd  of  2000  ppm  could  be  determined  for  its  mole  fraction. 

; 

/ 

The  flowrate  into  the  microprobe  varied  with  flame  position  due  to  changes  in  the  gas 
temperatmre  and  partial  clogging  of  the  orifice.  To  account  for  this,  all  PTMS  and  EQMS 
mass  spectra  were  normalized  with  the  EQMS  argon  signal  by  assuming  that  the  argon 
concentration  was  uniform  throughout  the  flame.  As  mentioned  above,  argon  was  added  to 
the  fuel  so  that  its  concentration  was  1%  in  both  reactant  streams.  At  a  height  close  to  the 
burner  siuface,  where  soot  clogging  was  insi^ificant,  the  radial  temperature  profile  could 
be  recovered  from  the  raw  argon  signals,  indicating  that  the  constant-argon  assumption  is 
valid. 

4.  Numerical  Method 

Our  goal  is  to  obtain  a  discrete  solution  of  the  two-dimensional,  elliptic,  governing 
conservation  equations  for  the  temperature,  gas-phase  and  soot  mass  fractions,  velocities 
and  vorticity  on  a  two-dimensional  mesh  the  initial  nodes  of  which  are  formed  by  the 
intersection  of  the  lines  of  the  mesh  Air 

Afr  =  {0  =  To  <  ri  <  . . .  <  r, . . .  <  r*/,  =  Rq},  (4.1) 

and  the  mesh  A4, 

Mz  =  {0  =  Zo  <  zi  <  . . .  <  Zj  . . .  <  zm.  =  2!}.  (4.2) 
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Computationally,  we  combine  a  steady-state  and  a  time-dependent  solution  method.  A 
time-dependent  approach  is  used  to  help  obtain  a  converged  numerical  solution  on  an 
initial  coarse  grid.  Grid  points  are  then  inserted  adaptively  and  Newton’s  method  is 
used  to  complete  the  problem.  Specifically,  the  spatial  operators  in  the  governing  partial 
differential  equations  are  approximated  by  finite  difference  expressions.  Diffusion  terms 
are  approximated  by  centered  differences  and  convective  terms  by  upwind  approximations. 
With  the  difference  equations  written  in  residual  form,  we  seek  the  solution  U*  of  the 
system  of  nonlinear  equations 

F(U)  =  0.  (4.3) 

For  an  initial  solution  estimate  U°  which  is  sufficiently  close  to  U*,  the  system  of  nonlinear 
equations  in  (4.3)  can  be  solved  by  Newton’s  method.  This  leads  to  the  iteration 

J(C7”)(17’'+' -  C/")  = -A”F(C/"),  n  =0,1,2,....  (4.4) 

where  J(U")  =  dF(U’')/dU  is  the  Jacobian  matrix  and  A"  (0  <  A  <  1)  is  a  damping 
parameter. 

We  point  out  that  with  the  spatial  discretizations  used  in  forming  (4.3),  the  Jacobian 
matrix  in  (4.4)  can  be  written  in  block-nine-diagona;l  form.  We  form  several  columns  of  the 
Jacobian  simultaneously  using  vector  function  evaluations  and  the  Jacobian’s  given  sparsity 
structure.  The  Newton  equations  are  solved  by  a  preconditioned  BiCGSTAB  iteration. 
The  Newton  iteration  continues  until  the  size  of  -  Z7"|l2  is  reduced  appropriately. 

Grid  points  of  the  two-dimensional  mesh  axe  determined  by  equidistributing  positive  weight 
functions  over  mesh  intervals  in  both  the  r  and  z  directions.  During  the  pseudo-transient 
portion  of  the  calculation,  the  size  of  the  time  steps  is  chosen  by  monitoring  the  local 
truncation  error  of  the  time  discretization  process.  Due  to  the  cost  of  forming  the  Jacobian 
matrices  with  detailed  transport  and  finite  rate  chemical  kinetics,  a  modified  Newton 
method  is  implemented  along  with  several  theoretical  estimates  that  determine  when  a  new 
Jacobian  should  be  reformed.  These  theoretical  results  help  increase  the  overall  efficiency 
of  the  algorithm  (see  Smooke  et  aJ.,  1990,  1996  for  more  details). 

5.  Results  and  Discussion 

In  this  section  we  apply  the  model,  solution  procedure,  and  experimental  methods 
discussed  in  Sections  2,  3  and  4  to  a  confined,  axisymmetric,  coflow,  methane-air  diffusion 
flame  with  the  goal  of  helping  to  interpret  flame  structure  and  soot  formation  phenom¬ 
ena.  Initial  modeling  results  were  compared  with  the  experimental  values  for  bulk  flame 
parameters;  i,e.,  flame  height,  peak  flame  temperatures,  etc.  Initially,  flame  height  was 
overpredicted  by  the  model  by  15%  and  temperatures  were  underpredicted  by  about  100- 
150K.  This  lack  of  agreement  was  considered  unacceptable.  Upon  analysis  of  both  the 
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model  and  the  experiment,  it  became  apparent  that  conduction/radiationfrom  the  base  of 
the  flame  to  the  burner  exit  tube  in  the  experiment  caused  noticeable  preheating  of  both 
inlet  fuel  and  air  streams.  Calculations  demonstrated  that  the  flame  height  and  local  tem¬ 
peratures  were  significantly  affected  by  these  assumed  inlet  temperatures.  Increasing  the 
inlet  temperatures  from  300  to  420K  (mass  flow  held  constant)  provided  better  agreement 
to  the  flame  height  and  local  temperatures  (see  later  discussion).  This  approach,  while 
not  perfect,  offers  a  tractable  computational  solution  to  a  challenging  problem.  This  in¬ 
crement  is  consistent  with  attempts  to  measure  the  temperature  at  the  burner  exit  plane. 
The  computational  results  demonstrate  the  importance  of  accurate  specification  of  ini¬ 
tial  conditions,  specifically  temperature  (and  resultant  effects  on  local  velocities),  on  the 
ntunerical  solutions. 

Time-stepping  from  nearly  converged  solutions,  led  to  the  build  up  of  large  soot  parti¬ 
cles  at  the  base  of  the  flame  in  the  wake  of  the  burner  lip.  Several  numerical  experiments 
demonstrated  that  the  source  of  this  build-up  was  thermophoresis.  Small  soot  particles 
generated  inside  the  flame  front  at  the  lip  were  driven  inside  this  small  wake  with  large 
temperature  gradients.  Neither  surface  growth  nor  inception  enhanced  the  soot  levels  in 
these  regions.  Coalescence  over  long  times  led  to  the  formation  of  large  particles  which  were 
not  flushed  away  by  the  flow  in  this  wake  region,  yet  were  continually  fed  as  time  elapsed 
by  small  particles  near  the  flame  front.  In  heavier  sooting,  co-flow  flames,  such  build-up 
with  time  has  been  observed  experimentally  as  evidenced  by  the  occasional  requirement 
to  shut  down  and  clean  particulates  from  the  exit  tube  (Santoro,  1996).  This  build-up  is 
not  observed  experimentally  for  the  lightly  sooting  methane  flames.  This  process  has  been 
suppressed'  in  the  model  by  setting  the  thermophoretic  velocity  to  zero  below  750K.  This 
does  not  cause  any  noticeable  change  of  predicted  soot  levels  in  the  remainder  of  the  flow 
field. 

The  chemical  kinetic  mechanism  used  in  this  work  was  derived  from  GRIMech  2.11 
with  its  recommended  thermodynamics  (Bowman  et  al.,  1995).  All  reactions  related  to 
NO  formation  and  other  odd  nitrogen  species  were  stripped  from  the  set.  In  addition,  a 
series  of  reactions  related  to  the  formation  and  oxidation  of  benzene,  and  related  species, 
were  included  (see  Table  1  of  Hall  et  aJ.,  1997).  The  dominate  ring  forming  reaction  was 
propargyl  recombination.  This  is  consistent  with  material  balance  analyses  of  benzene 
formation  mechanisms  based  on  the  experimental  species  concentrations.  The  burner  fuel 
tube  had  an  inner  radius  of  R/=0.5556  cm  with  a  wall  thickness  of  0.0794  cm.  The  radius 
of  the  coflow  was  Ro  —  4.7625  cm.  Conditions  at  the  fuel  and  oxidizer  jet  inlets  were  such 
that 

vf  =  5.52cm/sec,  vqx  =  12.54cm/sec,  (5.1) 

Tf  =  420K,  Tox  =  420K.  (5-2) 

All  calculations  were  performed  on  an  IBM  RS/6000  Model  590  computer.  In  the  com- 
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putations  to  be  presented  12  soot  size  classes  were  included  in  the  model.  Starting  from 
a  converged  solution  for  a  methane-air  flame  without  the  sectional  equations,  we  typi¬ 
cally  obtained  converged  solutions  for  the  complete  gas-soot  problem  in  several  hours  of 
computer  time. 

Temperature  Predictions 

Examination  of  the  curves  in  Figure  2  demonstrates  that  the  model  predicts  the  general 
shape  and  structure  of  the  flame  (location  of  pealc  flame  temperatures,  etc.)  reasonably 
well.  Temperature  is  perhaps  the  most  important  of  the  flame  parameters  as  temperature 
has  a  direct  influence  on  virtually  all  the  other  flame  properties  (e.g.,  flame  height  and 
width,  species  concentrations,  soot,  etc.)  A  more  careful  comparison  of  the  profiles  in 
Fig.  2  shows  that  the  model  underpredicts  the  experimental  measurements  along  the 
centerline  by  about  lOOK  at  all  heights  within  the  flame.  In  addition,  while  the  location 
and  value  of  the  peak  temperature  are  reasonably  well  described  up  to  about  a  height  of  z 
=  2  cm,  temperatures  are  overpredicted  at  the  outside  edge  (i.e.,  predicted  temperatmes 
do  not  fall  off  as  rapidly  as  the  data  with  increasing  radius).  Finally,  the  experimental 
flame  closes  to  the  centerline  (at  z  fa  40  mm)  more  rapidly  than  does  the  computed 
flame.  The  computed  low  centerline  temperatures  and  the  slightly  delayed  contraction 
of  the  flame  tip,  presumably  are  very  much  related  to  the  inability  to  match  accurately 
the  flame  height.  These  deficiencies  in  the  model  are  attributable,  at  least  partly,  to  the 
imcertainties  in  the  inlet  boundary  conditions.  As  argued  below,  these  limitations  of  the 
computed  temperature  field  have  a  direct  implication  on  the  prediction  of  gas  species  and 
particulates.  One  possible  explanation  of  the  low  temperatures  along  the  centerline  is  the 
lack  of  inclusion  in  the  model  of  absorption  by  methane  of  energy  radiated  from  the  flame 
front. 

Prediction  of  Species 

Figure  3  shows  comparison  of  the  methane  profiles  as  a  function  of  radius  at  several 
heights  above  the  burner.  Again  the  general  features  of  the  radial  profiles,  for  several 
heights  in  the  flame  are  reproduced  well  by  the  model.  Typically,  methane  is  overpredicted 
by  a  factor  of  two  along  the  centerline.  This  result  is  consistent  with  the  low  temperatures 
from  the  computed  flame  along  the  centerline.  Alternatively,  the  model  overpredicts  the 
fall-off  with  increasing  radius,  most  noticeably  at  the  lowest  height,  z  =  1.0  cm.  This 
accentuated  fall-off  is  undoubtedly  a  direct  result  of  the  higher  predicted  temperatures  in 
the  wings  of  the  flame,  and  hence  more  rapid  chemistry.  The  factor  of  three-overprediction 
at  z  =  2.5  cm  on  the  centerline,  is  a  consequence  of  the  (experimental)  rapid  decrease  in 
methane  concentration  in  this  region  of  the  flame,  i.e.,  the  flame  tip,  and  the  inability  of 
matching  the  flame  height  exactly. 

Similar  comments  can  be  made  in  regard  to  the  acetylene  and  benzene  profiles.  In 
these  cases,  however,  centerline  values  are  underpredicted,  presumably  because  methane 
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decomposition  is  suppressed  at  the  low  predicted  temperatures  along  the  centerline,  as 
discussed  above.  Not  surprisingly,  the  relatively  small  underprediction  of  acetylene  (factor 
of  two)  along  the  centerline,  leads  to  a  larger  deficit  in  the  computed  benzene  profiles  (a 
factor  of  three)  on  the  centerline.  Since  these  species  are  key  to  the  growth  and  incep¬ 
tion  of  particulates,  the  combination  of  their  low  values  as  well  as  the  low  temperature 
predictions  along  the  centerline  must  lead  commensurately  to  low  predictions  of  soot  in 
these  regions.  Of  importance,  however,  is  that  the  concentrations  of  these  species,  as  well 
as  temperatures,  are  better  described  in  the  flame  front  region  most  critical  to  inception 
and  growth.  Until  a  satisfactory  solution  is  found  to  predict  more  accurately  centerline 
temperatures,  quantitative  predictions  of  formation  of  the  translucent  droplets/particles 
along  the  centerline,  such  as  observed  by  Megaridis  and  Dobbins  (1989)  will  likely  not  be 
achieved. 

Soot  Particulates 

Figure  6  depicts  the  computed  and  experimental  values  of  soot  volmne  fraction  (cm^ 
soot/cm®  gas)  as  a  function  of  radius  at  2  =  1.75,  2.0.  2.25,  2.5,  2.75,  and  3.0  cm.  As  for 
the  species  profiles,  the  computations  describe  the  general  features  of  the  soot  field.  Peak 
soot  volume  fractions  are  underpredicted  by  a  factor  of  about  three,  however  the  shape 
and  trend  of  the  shift  of  the  volume  fraction  profile  is  reproduced  well.  The  experimental 
peak  values  collapse  to  the  centerline  more  rapidly,  as  discussed  above,  because  of  the 
difference  in  the  computed  fiame  height.  Centerline  values  are  significantly  underpredicted 
(order  of  magnitude)  because  of  the  combination  of  predicted  low  temperatures  (lOOK), 
low  acetylene  (factor  of  two)  and  low  benzene  (factor  of  three).  It  is  reasonable  to  assume 
that  should  these  bulk  flame  properties  be  more  accurately  modeled,  our  ability  to  attain 
quantitative  predictions  of  the  soot  field  would  be  enhanced  dramatically. 

Predicted  soot  voltune  fractions  are  strongly  dependent  on  predicted  acetylene  pro¬ 
files,  which  generally  are  dependent  upon  the  kinetics  supplied  by  GRIMech,  but  also  are' 
depressed  when  scrubbing  is  included  in  the  calculations.  Without  the  added  chemical 
mechanisms  and  soot  formation  steps,  predicted  levels  of  acetylene  are  very  close  to  the 
experimental  values,  depending  on  the  fiame  location.  Measured  peak  acetylene  concentra¬ 
tions  are  about  12,000  ppm.  The  measured  peak  soot  volume  fractions  are  4.7x10“^  which 
is  approximately  equivalent  to  4000  ppm  acetylene.  Hence,  suppression  of  the  acetylene 
peak  by  the  soot  formation  is  expected  and  is  observed  in  these  calculations  -  even  in 
these  lightly  sooting  flames.  The  reduced  acetylene  levels,  in  turn,  suppress  the  levels  of 
soot  predictions,  since  both  the  inception  and  the  growth  processes  are  strongly  depen¬ 
dent  on  the  acetylene  concentration.  Predicted  acetylene  concentrations  more  consistent 
with  experimental  results  would  result  in  more  favorable  agreement  between  predicted  and 
measured  soot  volume  fraction  in  Figure  6.  Soot  volume  fraction  predictions  based  on  this 
growth  model  have  some  residual  sensitivity  to  the  number  of  size  classes  at  the  value  of 
12  employed,  which  was  the  largest  number  permissible  given  the  number  of  species  and 
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grid  nodes  relative  to  available  computer  memory.  While  size  class  convergence  will  have 
to  be  addressed  ultimately,  factors  such  as  inlet/boundary  conditions  and  flame  chemistry 
may  have  a  much  greater  influence  on  the  accuracy  of  the  solution  and  its  ability  to  repro¬ 
duce  bulk  flame  parameters  (specifically,  temperature).  As  will  be  seen,  calculations  with 
another  soot  growth  model  show  much  more  rapid  convergence  with  the  number  of  size 
classes. 

2D  Solutions 

In  Figures  7-11  we  illustrate  contour  plots  of  the  same  variables  plotted  in  Figures  2-6. 
The  two-dimensional,  color  contour  of  computed  temperatures  in  Fig.  7  clearly  shows  the 
peak  temperature  regions  in  the  flame  front.  The  highest  temperatures  within  the  flame 
are  not  at  the  flame  peak,  but  in  the  outer  annulus  regions.  Peak  height  can  be  inferred 
from  these  profiles  by  locating  the  point  at  which  the  temperature  peaks  on  the  centerline, 
for  example.  The  contour  plot  of  methane  shows  a  very  dramatic  loss  of  methane  very  low 
in  the  flame.  This  extremely  rapid  decay  is  principally  due  to  rapid  diffusion  of  molecular 
nitrogen  to  the  centerline,  rather  than  the  reaction  of  methane.  (In  fact,  it  is  exactly  this 
rapid  diffusional  process  that  has  led  to  confusion  in  the  literature  when  attempts  have 
been  made  to  separate  dilution  from  temperature  effects  and  their  relative  importance  to 
soot  formation.).  Predicted  peak  acetylene  and  benzene  profiles  occur  on  the  centerline 
about  a  centimeter  below  the  flame  tip  near  z  =  2.8  cm.  Benzene  seems  to  decay  more 
rapidly  above  this  location  than  does  acetylene.  The  rapid  decay  in  the  acetylene  appears 
more  closely  aligned  with  the  flame  tip.  Alternatively,  the  earlier  decay  in  the  benzene 
is  hnked  tq  the  soot  inception  process,  which  attains  its  peak  rates  in  this  region  of  the 
flame. 

Soot.  Growth  and  Inception 

Figures  12-14  contain  contour  plots  of  the  net  soot  production,  soot  inception  and  soot 
oxidation,  respectively.  Comparison  of  the  net  soot  production  (Figure  12)  with  inception 
(Figure  13)  shows  that  the  growth  region  is  surface-  rather  than  inception-dominated. 
Higher  in  the  flame  net  growth  gives  way  to  oxidation  as  the  particles  burn  out  (Figure  14). 
The  surfa.ce  growth  contours  are  a  convolution  of  the  temperature  and  acetylene  contours 
(Figures  7  and  9),  strongly  reflecting  the  annular  shape  of  the  temperature  profile,  while 
inception  is  an  image  of  the  benzene  contour  of  Figure  10. 

Soot  Oxidation 

The  rates  of  oxidation  of  the  particulates  by  OH  and  by  O2  are  shown  in  Figure  14. 
Rate  data  are  plotted  on  the  same  scale  in  the  two  figures  to  readily  assess  the  relative 
importeince  of  the  two  oxidants.  Maximum  oxidation  rates  by  O2  clearly  are  about  a 
factor  of  5  to  10  below  those  by  OH.  Oxidation  by  O2  occurs  in  two  separate  regions  of 
the  flow.  Near  the  base  of  the  flame,  O2  diffuses  underneath  the  flame  and  contributes  to 
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soot  oxidation  on  the  fuel-rich  side  of  the  flame  until  the  O2  concentration  is  exhausted. 
O2  also  contributes  to  soot  oxidation  in  the  post  flame  region,  but  only  competes  with 
oxidation  by  OH  high  in  the  flame.  The  hydroxyl  radical  is  clearly  the  dominant  oxidizing 
species  in  this  flame.  Peah  oxidation  rates  occur  in  the  conventional  wing  tips  of  the 
flame  and  high  oxidation  rates  continue  to  the  centerline.  These  results  (i.e.,  the  strong 
dominance  of  OH  as  the  principal  oxidizing  species)  are  wholly  consistent  with  those  of 
Gaxo,  et  al.  (1990)  who  analyzed  the  dynamics  of  a  methane-air  coflow  diffusion  air  and 
more  recently  of  Kennedy,  et  al.  (1996)  in  their  detailed  modeling  of  an  ethylene  coflow 
diffusion  flame.  Together,  the  present  and  prior  results  cast  doubts  on  the  validity  of 
models  which  neglect  the  role  of  oxidation  by  OH  in  coflow  flames  (for  example,  Kaplan, 
et  al.  (1996)  and  Sivanthanu  and  Gore  (1994)).  The  cause  of  the  high  importance  of  OH 
as  the  principal  oxidizing  species  is  its  high  superequilibrium  level  in  the  flame  front.  (See 
discussion  below.) 

In  this  study,  the  soot  field  has  been  modeled  using  a  constant  0.13  collision  frequency 
for  reaction  after  Neoh,  et  al.  (1981).  In  the  Kennedy,  et  al.  (1996)  study,  this  reaction 
eflficiency  was  varied  linearly  from  0.05  low  in  the  flame  to  0.2  in  the  tip  flame  region, 
to  assist  in  matching  soot  profiles.  In  addition,  Garo,  et  al.  interpreted  their  data  to  be 
consistent  with  a  varying  collision  frequency  with  height  in  the  flame.  While  it  is  apparent 
that  such  methods  could  assist  matching  the  experimental  and  modeling  profiles  of  soot, 
we  have  not  yet  adopted  this  approach. 

Superequilibrium  OH 

Smook6,  et  al.  (1992)  have  computed  superequilibrium  levels  of  hydroxyl  in  methane 
air  coflow  flames.  Santoro  and  coworkers  have  measured  OH  concentrations  and  tempera¬ 
tures  in  ethylene  flames  and  found  these  radical  concentrations  to  be  effected  by  the  level 
of  soot  in  the  flame.  Higher  soot  levels  decreased  OH  concentrations  due  to  scrubbing. 
Such  effects  were  not  noticeable  in  this  study,  presumably  due  to  the  much  lower  concen¬ 
trations  of  soot.  The  OH  measured  by  Puri,  Semtoro  and  Smyth  (1994)  were  described 
to  be  at  superequilibrium  levels.  Logically,  these  superequilibrium  levels  are  the  source  of 
the  artificial  requirement  to  augment  the  reaction  efficiency  of  OH  and  soot.  To  exam¬ 
ine  levels  of  superequilibrium,  calculated  OH  (referred  to  here  as  “flame  OH”  have  been 
compared  to  the  equiUbrium  levels  of  OH  based  on  1)  the  local  temperature  (constant 
temperature  calculation),  2)  the  full  equilibrium  (constant  enthalpy  calculation),  and  3)  a 
partial  equilibrium  assumption  based  on 

02  +  Hi<-*20H 

and  the  local  flame  temperature.  The  ratio  of  this  flame  OH  to  the  equilibrium  levels 
are  plotted  in  Figure  15  as  a  function  of  the  radius  at  ^  =  4.01  cm.  The  maximum  soot 
oxidation  rate  occurs  at  this  height  at  a  radius  of  about  0.31  cm.  At  this  radius,  the 
OH  attains  concentrations  a  faetor  about  10  above  the  equilibrium  value,  thus  validating 
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further  the  existence  of  large  OH  superequilibrium  concentrations  as  an  important  con¬ 
tributor  to  the  destruction/burnout  of  soot.  Chemical  systems  which  cannot  support  such 
high  superequilibrium  levels  may  be  more  prone  to  soot  emissions.  The  partial  equilibrium 
assumption  (assuming  knowledge  of  the  H2,  02  concentrations  and  temperature)  provides 
a  much  better  estimate  of  hydroxyl,  at  least  through  the  flame  front  region  where  soot  is 
being  oxidized. 

Partcle  Size 

Figure  16  shows  the  evolution  of  soot  particle  size  distribution  in  space  as  a  result  of 
surface  growth  and  coalescence.  The  sectional  mass  boundaries  have  been  placed  linearly 
on  a  logarithmic  scale.  Peak;  soot  mass  concentration  occurs  in  size  class  nine,  whose 
mass-weighted  mean  diameter  is  about  48  nm,  which  should  be  considered  an  estimate.  A 
larger  number  of  sections  than  that  employed  here  would  be  required  to  obtain  an  accurate 
average  particle  size.  This  average  size  is  clearly  larger  than  typical  determinations  of  20- 
30  nm  in  laboratory  flames  and  suggests  some  refinements  in  our  soot  model  are  still 
required.  The  fact  that  the  suspended  mass  in  the  highest  class  remains  relatively  small 
shows  that  the  upper  boundary  of  the  particle  size  domain  has  been  set  high  enough  to 
avoid  unphysical  accumulation  of  particles  in  the  highest  size  class.  The  smallest  size 
classes  are  observed  to  disappear  rapidly  in  the  post  flame  zone,  while  a  portion  of  the 
largest  particles  survive  to  the  end  of  the  computational  domain. 

Particle  size  distribution  within  the  flame  can  also  be  visualized  from  Figures  17  and 
18  in  which  the  mass  fractions  of  each  of  the  12  soot  mass  classes  (me)  are  plotted  as  a 
function  of  radial  distance  for  heights  of  1  and  3  centimeters.  Low  in  the  flame,  the  sixth 
soot  class  contains  the  peak  mass  fraction  and  at  3  cm  (just  before  the  onset  of  rapid  soot 
oidation  by  OH  inside  the  flame  front),  the  maximum  occurs  in  the  ninth  class.  At  this 
higher  position,  substantial  soot  is  found  also  near  the  centerline  of  the  flame.  However, 
the  average  particle  for  this  “centerline”  soot  is  much  smaller  than  that  just  inside  the 
flame  front,  a  result  qualitatively  consistent  with  the  results  by  Megaridis  and  Dobbins 
(1989).  Since  benzene  concentrations  (Figure  10)  and  hence  inception  rates  (see  Figure 
13)  are  peaking  in  this  region  of  the  flow,  the  particles  depicted  on  the  centerline  are  likely 
formed  along  or  near  to  the  centerline  rather  than  transported  via  thermophoresis  from 
the  flame  front  to  this  centerline  region. 

Radiation  Effects 

While  it  might  be  expected  that  thermal  radiation  effects  would  be  sinall  in  a  lightly 
sooting  flame  like  that  investigated  here,  simulations  in  which  the  radiative  loss  term,  Eq. 
2.21,  is  set  equal  to  zero  show  that  radiation  effects  are  significant.  Without  radiation, 
predicted  peak  temperatures  are  122K  higher,  a  change  large  enough  to  affect  predicted 
soot  levels  and  species  concentrations.  Predicted  peak  volume  fractions  increase  by  a 
factor  of  three,  with  peak  acetylene  and  benezene  levels  increasing  by  45%  and  22%, 
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respectively.  The  significance  of  thermal  radiation  in  this  lightly  sooting  fiame  arises  from 
residence  times  that  are  relatively  long  compared  to  those  for  counterfiow  diffusion  flames, 
in  which  radiation  effects  are  significant  only  for  extremely  low  strain  rates  (Garcia,  et 
al.,  1996).  With  the  relatively  low  values  of  soot  concentration  and  the  fact  that  the 
soot  does  not  occupy  the  peak  temperature  regions,  the  thermal  radiation  here  is  actually 
dominated  by  water  vapor  and  carbon  dioxide,  with  the  4.3  micron  band  of  CO2  the 
dominant  contributor.  The  physical  dimensions  of  this  flame  are  sufficiently  small  that 
self-absorption  or  optical  thickness  corrections  should  be  minimal. 

Variations  of  Inlet  Gas  Temperatures 

Early  modeling  efforts  with  the  inlet  gas  temperatures  at  298K  resulted  in  a  flame 
height  about  15%  higher  than  the  experimental  value.  In  addition,  predicted  peak  flame 
temperatures  were  about  125K  below  the  experimental  value.  Since  soot  volume  fractions 
axe  a  strong  function  of  temperatures,  these  low  predicted  temperatures  substantially  lim¬ 
ited  our  ability  to  predict  accurately  the  soot  concentrations.  Local  temperatures  in  the 
unburned  gases  just  above  the  burner  lip  were  measured  to  be  much  higher  than  298K, 
presumably  due  to  radiation  preheating  and  conduction.  Accurate  treatment  of  this  prob¬ 
lem  would  require  knowledge  not  only  of  temperatures,  but  also  of  local  velocities  at  each 
location  in  the  burner  inlet.  To  1)  analyze  this  problem  in  a  simplified  manner,  2)  provide 
information  relative  to  a  sensitivity  study,  and  3)  elevate  the  peak  flame  temperatures 
to  values  closer  to  the  experimental  levels,  calculations  were  performed  for  elevated  tem¬ 
peratures  in  the  entire  fuel  flow,  the  entire  air  flow,  or  both  streams.  For  each  case,  the 
mass  flowS;  per  unit  area  were  held  constant  and  hence  velocities  increased  with  higher 
temperatures. 

Table  1  compares  the  peak  soot  volume  fractions  and  the  peak  flame  temperatures 
obtained  for  several  solutions  at  inlet  temperatures  ranging  from  298  to  500K.  An  ex¬ 
amination  of  these  numbers  demonstrates  that  an  increase  in  the  air  temperature  has  a 
much  more  dramatic  effect  on  both  the  peak  flame  temperature  and  volume  fraction  than 
comparable  changes  in  the  fuel  temperature.  Presumably,  this  effect  is  due  to  the  higher 
mass  weighted  enthalpy  content  of  air  under  stoichiometric  conditions.  The  higher  flame 
temperatures  which  arise  from  the  elevated  inlet  air  temperatures  lead  to  a  proportionally 
greater  increase  in  the  soot  generation. 

Comparisons  to  Fairweather,  et  al.  model 

Finally,  the  model  was  rerun  with  the  inception  and  surface  growth  model  of  Fair- 
weather  et  al.  (1992).  Soot  volume  fraction,  net  soot  production,  inception,  and  oxidation 
are  illustrated  in  Figures  19-22.  In  this  model,  which  has  been  applied  to  ethylene-  and 
acetylene^fuelled  coflow  flames  (Kennedy,  1996;  Sivathanu  and  Gore,  1994),  both  inception 
and  surface  growth  are  proportional  to  acetylene  concentration.  Its  inception  and  surface 
growth  rates  are  used  here  in  the  polydisperse,  sectional  representation  rather  than  the 


19 


monodisperse  representation  of  the  original  model.  Additionally,  oxidation  by  OH,  not 
included  in  Fairweather,  et  al.  (1992),  has  been  accounted  for  in  this  implementation. 


TABLE  1 

Effect  of  Fuel  and  Oxidizer  Inlet  Temperatures  on 
Peak  Soot  Volume  Fraction  and  Peak  Flame  Temperature 


Tj/Z, 

298/298 

3.30  X  10-® 

1942  K 

298/420 

8.71  X  10-® 

1976  K 

298/500 

1.43  X  10"^ 

2000K 

420/298 

4.67  X  10"® 

1948  K 

500/298 

5.87  X  10-* 

1954  K 

420/420 

1.14  X  10-"^ 

1984  K 

500/500 

x.xx  X  lO"’^ 

XX  K 

The  predicted  soot  concentrations  (Figure  19)  are  seen  to  be  slightly  larger  than  those 
predicted  with  the  base  model.  Given  the  previously  discussed  uncertainties  associated 
with  the  chemical  mechanism  and  burner  boundary  conditions,  however,  no  judgements 
about  the  merits  of  this  model  vis-a-vis  the  base  model  are  warranted.  Inception  is  more 
significant  in  the  Fairweather,  et  al.  (1992)  model  (compare  Figures  21  and  13),  and  this 
seems  to  give  calculations  using  it  much  less  sensitivity  to  the  number  of  soot  sections 
employed.  Predicted  soot  volume  fractions  using  five  and  twelve  sections  were  not  signifi¬ 
cantly  different.  The  spatial  distribution  of  inception  in  the  two  models  is  very  different, 
as  Figures  13  and  21  show.  The  Fairweather,  et  al.  (1992)  model  has  an  inception  profile 
that  is  strongly  annular,  reflecting,  like  the  surface  growth,  the  convolution  of  temperature 
and  acetylene  profiles.  Oxidation  in  this  latter  model  is  also  dominated  by  OH,  as  Figure 
22  illustrates. 

Summary  of  Differences  Betweeen  Model  and  Experiment 

Differences  between  our  predictions  and  the  experimental  data  include:  overprediction 
of  flame  height,  underprediction  of  peak  temperatures,  too  slow  decay  of  methane  along 
the  centerline  and  a  slow  decay  of  temperature  in  the  wings  of  the  flame.  In  addition, 
the  model  significantly  underpredicts  the  peak  soot  volume  fractions  in  the  flame  and 
tends  to  overpredict  particle  size.  Much  of  the  differences  can  be  directly  linked  to  the 
inadequate  prediction  of  the  temperature  field  and  the  flame  height.  Contrasting  these 
results  are  the  very  good  agreement  of  the  temperature  field  and  flame  height  obtained  in 
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the  modeling  of  a  methane  coflow  flame  in  a  study  of  NO  formation  Smooke,  et  al.  (1997). 
Principal  differences  between  this  previous  study  and  that  of  the  present  investigation 
includes  dilution  of  the  fuel,  higher  flow  rates,  and  a  lifted  flame.  We  believe  the  most 
likely  of  these  parameters  affecting  uncertainties  in  the  present  study  is  the  inability  of 
our  code  to  properly  model  the  attachment  (i.e.,  close  proximity)  of  the  flame  to  the 
burner  lip  and  all  of  its  attendant  effects.  Other  possible  contributors  to  this  error  include 
uncertainties  in  the  kinetics.  For  example,  using  this  kinetic  mechanism  (hydrocarbon 
portion  of  GRIMech  2.11)  the  extinction  limit  in  a  counterflow  methane-air  diffusion  flame 
is  significantly  overpredicted  (380  cm“^  experiments  vs.  predictions  of  620  cm“^).  We  note 
also  that  there  are  still  fundamental  uncertainties  in  soot  growth  modeling  such  as  the  rates 
of  inception,  surface  growth,  and  coalescence  of  very  small  particles.  For  example,  a  50% 
increase  in  the  surface  growth  rate  constant  results  in  a  predicted  soot  volume  fraction  in 
much  better  agreement  with  data.  Until  the  temperature  field  is  more  accurately  described, 
it  is  difficult  to  offer  quantitative  arguements  about  the  accuracy  (or  lack  thereof)  of  our 
soot  model.  Nevertheless,  the  qualitative  results  and  conclusions  should  still  be  valid. 

6.  Conclusions 

The  axisymmetric  laminar  diffusion  flame  provides  a  natural  environment  in  which  one 
can  investigate  the  interaction  of  soot  formation  with  detailed  gas-phase  chemistry  in  a 
multidimensional  system.  In  this  paper  we  have  compared  experimental  measurements  of 
the  temperature,  soot  volume  fraction  and  selected  species  from  an  axisymmetric,  lami¬ 
nar,  methetoe-air  diffusion  flame  with  computational  results  obtained  by  generalizing  the 
sectional  model  presented  in  Hall  et  al.  (1997)  to  a  coflow  model.  This  work  represents  the 
first  time  that  a  detailed  chemistry /complex  transport  flame  model  has  been  coupled  to  a 
polydisperse  soot  growth  analysis  in  which  the  scrubbing  of  gaseous  species  and  both  gas 
and  particle  radiation  are  included.  The  detailed  chemistry  includes  species  up  through 
benzene  and  includes  a  steady-state  model  for  the  formation  of  polyaromatic  hydrocarbons 
to  simulate  soot  inception.  Soot  volume  fractions  calculated  from  this  model  using  either 
our  growth  inception  formulation  or  that  by  Fairweather  et  al.  (1992)  underpredicted  the 
experimental  values  by  about  a  factor  of  three.  These  low  values,  we  believe,  are  directly 
attributable  to  the  inability  to  reproduce  accurately  bulk  flame  properties.  Specifically, 
low  temperatures  lead  to  underpredictions  of  critical  growth  species  such  as  acetylene  and 
benzene.  Principal  differences  between  the  model  and  the  experiments,  .we  believe,  are 
attributable  to  the  inaccuracies  in  the  predicted  temperature  field  which,  in  turn,  may  be 
partially  due  to  inadequate  modeling  of  flame  attachment  at  the  burner  lip. 
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Figure  Captions 


Figure  1  Schematic  of  the  burner  configuration. 

Figure  2  Profiles  of  the  computed  and  measured  temperature  profiles  at  various  heights 
in  the  flame. 

Figure  3  Profiles  of  the  computed  and  measured  methane  profiles  at  various  heights  in 
the  flame. 

Figure  4  Profiles  of  the  computed  and  measured  acetylene  profiles  at  various  heights  in 
the  flame. 

Figure  5  Profiles  of  the  computed  and  measured  benzene  profiles  at  various  heights  in 
the  flame. 

Figure  6  Profiles  of  the  computed  and  measured  soot  volume  fraction  profiles  at  various 
heights  in  the  flame. 

Figure  7  Temperature  isotherms  of  the  computational  model  as  a  function  of  the  radial 
and  axial  coordinates. 

Figure  8  Methane  isopleths  of  the  computational  model  as  a  function  of  the  radial  and 
axial  coordinates. 

Figure  9  .^cetylene  isopleths  of  the  computational  model  as  a  function  of  the  radial  and 
axial  coordinates. 

Figure  10  Benzene  isopleths  of  the  computational  model  as  a  function  of  the  radial  and 
axial  coordinates. 

Figure  11  Soot  volume  fraction  isopleths  of  the  computational  model  as  a  function  of  the 
radial  and  axial  coordinates. 

Figure  12  Net  soot  production  isopleths  of  the  computational  model  as  a  function  of  the 
radial  and  axial  coordinates. 

Figure  13  Comparison  of  soot  inception  due  to  acetylene  and  benzene  as  a  function  of 
the  radial  and  axial  coordinates;  rate  of  equation  2.1  versus  that  of  equation  2.2. 

Figure  14  Comparison  of  soot  oxidation  due  to  oxygen  and  OH  eis  a  function  of  the  radial 
and  axial  coordinates. 

Figure  15  Profiles  of  the  ratio  of  the  OH  mole  fraction  to  the  equilibrium  levels  as  a 
function  of  the  radial  coordinate  at  z  =  4.01  cm. 
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Figure  16  Illustration  of  the  evolution  of  soot  particle  size  distribution  in  space  as  a 
result  of  surface  growth  and  coalescence.  The  sectional  mass  boundaries  have  been  placed 
linearly  on  a  logarithmic  scale.  Peak  soot  mass  concentration  occurs  in  size  class  nine, 
whose  mass- weighted  mean  diameter  is  about  48  nm. 

Figure  17  Contributions  of  the  Soot  Mass  Classes  to  the  Total  Soot  Volume  Fraction  at 
z=l  cm. 

Figure  18  Contributions  of  the  Soot  Mass  Classes  to  the  Total  Soot  Volume  Fraction  at 
z=3  cm. 

Figure  19  Soot  volume  fraction  isopleths  of  the  Fairweather  et  al.  (1992)  model  as  a 
function  of  the  radial  and  axial  coordinates. 

Figure  20  Net  soot  production  isopleths  of  the  Fairweather  et  al.  (1992)  model  as  a 
function  of  the  radial  and  axial  coordinates. 

Figure  21  Soot  inception  isopleths  of  the  Fairweather  et  al.  (1992)  model  as  a  function 
of  the  radial  and  axial  coordinates. 

Figure  22  Comparison  of  soot  oxidation  due  to  oxygen  and  OH  as  a  function  of  the  radial 
and  axial  coordinates  for  the  Fairweather  et  al.  (1992)  model. 
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Interpretations  of  a  Computational  Study  of  Soot  Formation 
a  Coflow,  Laminar  Diffusion  Flame 
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Introduction 

Coflow  axisymmetric  flames  are  often  employed  as  a  base  flame  for  describing  practical 
combustion  devices  as  well  as  for  gaining  information  on  interactions  between  heat  and  mass 
transfer  as  well  as  detailed  chemistry.  The  ability  to  predict  these  coupled  effects  is  critical  to  the 
modeling  of  turbulent  reacting  flows  and  to  the  understanding  of  the  controlling  phenomena 
related  to  generation  of  pollutants.  Most  studies  using  detailed  chemical  kinetics  and  coupled 
models  of  soot  production  and  oxidation  for  examining  computational  combustion  have  focused 
on  one  dimensional  geometries.  A  few  recent  studies  [e.g.,  1,2]  have  focused  on  multi¬ 
dimensional  problems  in  which  monodisperse  soot  formation  models  and  skeletal  chemical 
kinetic  mechanisms  were  used.  This  abstract  describes  the  results  [3]  of  a  more  recent  modeling 
study  which  includes  both  detailed  chemical  kinetics  and  a  sectional  analysis  of  particle  dynamics 
and  then  presents  some  interpretations  of  model  predictions  for  a  lightly  sooting,  coflow  methane 
flame. 

Brief  Description  of  Model  and  Results 

In  the  recent  study  [3],  a  detailed  model  with  folly  coupled,  elliptic,  two-dimensional  equations 
including  soot  production  was  developed  and  compared  to  experimental  data  from  an 
axisymmetric,  laminar,  methane/air  coflow  diffusion  flame.  The  two-dimensional  model  includes 
a  detailed  transport,  finite  rate  (detailed)  chemistry  gas  phase  system  coupled  with  aerosol 
dynamical  equations  in  the  sectional  representation.  Detailed  submodels  for  the  inception,  growth, 
oxidation  and  dynamics  of  soot  particulates  as  well  as  effects  of  radiation  were  included  in  the 
model  [4,5].  Predictions  and  measurements  of  temperature,  soot  volume  fractions  and  selected 
species  were  compared  over  a  range  of  heights  and  as  a  function  of  radius.  Consistent  with  the 
experimental  data,  peak  soot  volume  fractions  were  predicted  in  an  annulus  region  just  inside  the 
flame  front.  If  radiation  losses  were  not  included,  the  peak  flame  temperature  increased  by  more 
than  lOOK  with  an  attendant  increase  by  a  factor  of  three  in  the  predicted  peak  soot  volume 
fractions.  Initial  modeling  results  significantly  underpredicted  the  peak  temperatures  and  soot 
volume  fractions  and  overpredicted  the  flame  height.  Increasing  gas  inlet  temperatures, 
particularly  the  air  temperature,  led  to  better  agreement  with  the  experimental  data  and  a  constant 


inlet  temperature  of  420K  was  used  for  additional  analysis.  The  model  predicted  the  flame  to  be 
about  4.8  cm  tall  (as  defined  by  peak  temperature  along  the  centerline).  Soot  growth  rates  start  to 
become  large  in  the  aimulus  region  of  the  flame  at  a  height  of  about  0.9  cm.  Peak  growth  rates  are 
found  at  about  2.2  cm  high  and  soot  undergoes  significant  decay  with  rapid  oxidation  occurring 
above  3.7  cm.  Many  interpretations  of  the  model  and  comparisons  to  data  are  presented  in  [3]. 

Interpretations 

Hydroxyl  radicals  in  the  flame  front,  and  most  importantly  in  the  soot  oxidation  region,  attain 
levels  which  are  about  10  times  higher  than  the  local  equilibrium  conditions.  A  comparison  of  the 
predicted  [OH]  as  a  function  of  radius,  relative  to  various  methods  for  estimating  the  radical 
concentration,  is  shown  in  Fig.  1  at  a  height  of  4  cm.  At  this  height  and  at  a  radius  of  0.3  cm,  peak 
soot  volume  fractions  are  rapidly  decaying  as  soot  rapidly  oxidizes  due  to  high  reaction  rates  with 
hydroxyl  radicals.  These  results  strongly  support  the  dominant  role  of  hydroxyl  radical  to  the 
oxidation  of  soot  in  coflow  flames  and  demonstrate  the  importance  of  super-equilibrium  [OH]. 
The  erroneous  assumption  of  equilibrium  [OH]  can  readily  lead  one  to  conclude  incorrectly  that 
treatment  of  oxidation  by  Oj,  alone,  is  sufficient. 

The  environment  in  which  soot  is  formed  was  examined  by  plotting  in  Fig.  2  the  locii  of  peak  soot 
surface  growth  rates  (at  each  height)  on  curves  of  temperature  as  a  function  of  radius.  Soot 
formation  (i.e.,  surface  grovlh  rates)  maximizes  on  the  rich  side  of  the  flame  at  temperatures  of 
1720-1760K,  a  temperature  which  is  in  the  vicinity  of  the  traditional  peak  of  the  well-known,  bell¬ 
shaped  curve  for  soot  formation.  The  lower  end  of  this  temperature  range  is  probably  a  better 
characteristic  temperature  for  the  peak  surface  growth  rate,  as  the  1760K  value  at  the  centerline, 
near  the  top  of  the  flame,  is  the  lowest  temperature  on  the  rich  side  attained  at  this  height  in  the 
flame.  Additional  information  about  the  soot  formation  region  can  be  extracted  firom  calculations 
of  the  local  equivalence  ratio,  based  the  local  gas  phase  C,  H,  O  elemental  concentrations.  These 
local  equivalence  ratios  are  plotted  in  Fig.  3,  as  a  function  of  radius  and  for  several  flame  heights. 
The  locii  of  peak  flame  temperatures  are  shown  and  confirms  that  peak  temperatures  occur  at  or 
near  the  rich  side  of  the  stoichiometric  flame  surface.  The  locations  of  peak  surface  groAvih  rates, 
as  well  as  those  of  peak  soot  volume  firactions,  are  also  plotted  in  Fig.  3.  Lower  in  the  flame,  these 
regions  coincide  near  equivalence  ratios  of  1.8,  separate  and  move  to  less  rich  regions  of  the  flow 
with  increasing  flame  height.  The  coincidence  of  the  peak  surface  growth  rates  and  the  peak  soot 
volume  fractions  is  a  significantly  different  result  from  that  of  counterflow  flames.  In  the  latter 
flames,  small  soot  particles  formed  at  high  temperatures  near  the  flame  front  are  rapidly  pushed 
[5]  by  both  thermophoresis  and  convection  towards  the  cooler  stagnation  plane  (about  900K) 
where  coagulation  occurs  readily  and  surface  growth  is  slow  due  to  the  low  local  temperatures. 
Thus,  it  is  likely  that  there  may  be  distinct  issues,  particularly  with  regard  to  the  relative 
importance  of  surface  growth  and  coagulation/coalescence  between  the  two  flame  types.  As  an 
example,  ‘ageing’  effects  should  be  more  significant  for  coflow  flames  in  which  soot  particulates, 
once  formed,  remain  in  a  high  temperature  zone.  Although  inception,  surface  growth  and 
coalescence  are  treated  in  our  model,  neither  ageing  nor  agglomeration,  however,  are  yet  included 
in  our  analysis  and,  hence,  the  impact  of  these  processes  on  model  solutions  is  uncertain. 

Also  of  interest  in  this  analysis  is  a  calculation  showing  the  effects  of  differential  species 
diffusion.  Radial  variations  of  the  local  H/C  molar  ratio  are  plotted  in  Fig.  4  for  several  flame 


heights.  These  calculations  clearly  show  the  effect  of  light  molecular  weight  species,  i.e.,  H  and 
H2  (and  presumably  some  HjO),  that  rapidly  diffuse  through  the  flame  front  and  out  beyond  the 
wings  of  the  flame.  A  resultant  minimum  in  the  H/C  ratio  generally  occurs  just  inside  the  peak 
flame  temperature  region  (shown  in  the  figure)  and  then  it  grows  again  with  increasing  radius  as 
there  is  no  carbon  in  the  coflow  air  stream.  The  impact  on  the  soot  formation  process  is  that  the 
local  H/C  ratios  where  soot  is  produced  are  noticeably  lower  (by  about  10%)  from  that  of  the  neat 
fuel. 

In  these  calculations,  12  size  classes  were  employed,  with  lower  and  upper  mass  limits  of  127  to 
1.5x10'®  amu.  This  upper  mass  boundary  (not  average  mass)  of  the  largest  size  class  was  picked 
such  that  the  solution  was  independent  of  this  boundary.  Particle  size  distributions  can  be 
determined  at  each  location  within  the  flame  and  are  not  a  priori  assumed.  Peak  volume  fraction 
occurs  in  size  class  nine,  with  a  mass-weighted  mean  diameter  of  48  nm.  This  relatively  large  size 
is  likely  larger  than  the  experimental  value,  perhaps  in  part  due  to  the  lack  of  ageing  effects  in  our 
model  formulation.  The  average  size  of  particles  is  much  lower  along  the  centerline,  consistent 
with  the  recent  studies  by  Dobbins  and  coworkers  [6]  in  an  ethylene  diffusion  flame.  This  flame 
would  normally  be  classified  as  non-sooting  in  that  the  experiment  indicates  that  soot  does  not 
escape  from  the  ‘wing-tips’.  A  close  examination  of  the  computational  solution  indicates  that 
nearly  all  the  soot  particles  are  oxidized  as  the  soot  field  passes  through  the  flame  front;  however, 
while  medium  and  smaller  particles  are  oxidized  readily  due  to  their  larger  surface-to-volume 
ratio,  small  concentrations  of  the  largest  particles  appear  to  escape  and  travel  downstream.  Our 
model  likely  overestimates  this  effect,  again  due  to  its  lack  of  inclusion  of  ageing  phenomena;  yet, 
even  so,  the  predicted  levels  are  likely  undetectable  as  they  are  about  40  times  lower  than  the  peak 
volume  fraction  elsewhere  in  this  lightly  sooting  flame. 
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